UNIVERSITY  OF  BASEL 


FRS  Klingelbergstrasse  80  Tel. 

Fax  +41  61  267  38  55 
E-Mail  J. P .Maier@unibas . ch 

http :  /  /www .  chemie  .  unibas  .  ch/~maier 


Department  of  Chemistry 

+  41 


Prof.  Dr.  John  P.  Maier, 
61  267  38  26  (or  00) 

CH-4056  Basel 
Switzerland 


Final  report:  IRI  Project:  contract  no.  F61775-02-C4029,  June  2002-May  2005 

The  broad  objectives  of  this  IRI  project  which  ran  during  three  years,  June  2002-  May  2005,  were  the 
investigations  of  the  spectroscopic  properties  of  boron  species,  their  aggregation  and  reactions  in 
matrices  ranging  from  neon  to  those  doped  with  hydrogen.  The  studies  of  the  aggregation  properties 
were  hampered  by  the  lack  of  spectroscopic  knowledge  on  the  electronic  transitions  of  the  polyatomic 
boron  molecules  and  their  ions.  Thus  it  was  clear  that  their  identification  at  first  had  to  be  undertaken 
in  order  to  proceed  on  this.  The  same  applies  to  reactions  in  hydrogen  dopped  matrices.  This  was 
accomplished  for  B3  and  B3' .  In  addition  theoretical  predictions  for  the  future  studies  of  the  B4  species 
and  its  positive  and  negative  ions  were  completed.  The  project  was  initially  instigated  on  the  advice 
and  recommendation  of  Dr.C.W.Larson,  AFRL,  Edwards  Air  Force  Base,  California,  whose  continual 
guidance,  suggestions  and  enthusiasm  throughout  the  project  were  much  appreciated. 

For  the  purpose  of  the  spectral  identification  and  characterisation  of  polyatomic  boron  molecules,  a 
world-wide  unique  apparatus,  combining  mass-selection  (using  a  cesium  sputter  source)  with  a  matrix 
isolation  technique  was  initially  used,  backed  up  by  theoretical  calculations.  The  six,  6-month  interim 
reports  summarized  the  experiments  undertaken,  including  the  unsuccessful  ones.  Because  the 
electronic  spectra  of  the  triatomic  boron  species  were  identified  in  5  K  neon  matrices,  an  ambitious 
project  to  observe  the  gas  phase  spectrum  of  B3  could  also  be  realized. 

The  project  as  a  whole  proved  rather  successful  and  productive  as  evidenced  by  5  publications  during 
this  three  year  period  in  leading  journals  of  chemical  physics: 

“Electronic  absoiption  spectra  of  B3  and  B3'  in  neon  matrices  and  ab  initio  analysis  of  the  vibronic 
structure” 

M.Wyss,  E.Riaplov,  A.Batalov,  J.P.Maier,  T.Weber,  W.Meyer  and  P.Rosmus 
Journal  of  Chemical  Physics,  119,  9703,  2003. 

“Gas  phase  detection  of  cyclic  B3 :  2  E  -  X  2  A3  electronic  origin  band“ 

P.Cias,  M.Araki,  A.Denisov  and  J.P.Maier 
Journal  of  Chemical  Physics,  121,  6776,  2004. 
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“The  near  infrared  l2  A2  -X  2Ai  electronic  transition  of  B3in  a  neon  matrix” 

A.Batalov,  J.Fulara,  I.Shnitko  and  J.P.Maier 
Chemical  Physics  Letters,  404,  315,  2005. 

“Vibrations  in  the  B4  rhombic  structure” 

R.  Linguerri,  I.Navizet,  P.Rosmus,  S. Carter  and  J.P.Maier 
Journal  of  Chemical  Physics,  122,  34301,  2005. 

“Theoretical  study  of  the  electronic  excited  states  of  B4 ,  B4+  and  Bf  “ 

C.Gillery,  R.Linguerri,  P.Rosmus  and  J.P.Maier 
Z.  Phys.Chem.,  219,  467,  2005. 

These  five  publications  summarize  the  results  and  success  of  this  IRI  project  and  are  included  with 
this  report. 

The  project  involved  experimental  studies  using  our  special  matrix  apparatus  in  Basel  with  which  a 
number  of  PhD  students  worked  during  this  period,  starting  with  M.Wyss,  E.Riaplov  (both  PhD  now) 
and  continued  by  A.Batalov,  I.Shnitko.  They  were  assisted  in  the  experiments  by  the  annual  visits  of 
Dr.  J.Fulara,  a  long  standing  collaborator  from  the  Polish  Academy  of  Sciences,  Warsaw,  and 
succeeded  in  the  measurement  of  the  electronic  spectra  of  B3  and  B3'  in  neon  matrices.  As  a  follow  up 
to  this  the  gas  phase  spectrum  of  B3  was  obtained  by  the  students  P.Cias  (now  PhD),  A.Denisov, 
working  under  the  guidance  of  Dr.M.  Araki. 

The  interpretation  of  the  vibronic  structure  and  symmetry  of  the  electronic  transitions  of  B3  and  B  3' 
proved  possible  as  result  of  extensive  collaborations  with  two  theoretical  groups:  that  of 
Prof.P.Rosmus,  University  of  Marne  la  Vallee,  France  and  Prof.W.Meyer,  University  of 
Kaiserslautern,  Germany.  P.Rosmus  spent  1-2  month  periods  in  each  of  the  three  years  of  this  IRI 
project  as  a  guest  Professor  here  working  on  the  structure  and  spectroscopy  of  boron  species, 
predicting  their  geometric,  electronic  and  vibrational  patterns.  Two  of  his  PhD  students  (C. Gillary, 
R.Linguerri)  also  had  sojourns  in  Basel  and  these  theoretical  studies  are  a  major  part  of  their  theses. 

Though  the  IRI  project  has  come  to  an  end,  we  plan  to  continue  our  spectroscopic  studies  of  the  boron 
clusters.  Specifically  the  tetraatomic  and  larger  boron  molecules  could  be  studied  in  the  gas  phase  by 
multiphoton  ionisation  and  laser  induced  fluorescence  methods-  the  instrumentation  for  which  is 
available  in  our  laboratory  in  Basel-  combined  with  a  laser  ablation  source.  For  the  neon  matrix 
experiments  a  new  ion  source,  incoiporating  a  high  repetition  laser  vaporization  design,  for  the 
deposition  of  the  larger  boron  ions  is  in  construction. 

EOARD  is  thanked  gratefully  for  providing  this  three  year  contract  which  has  enabled  us  to  open  up  a 
new  area  of  research  for  boron  clusters,  led  to  multi-national  collaborations  and  resulted  in  interactions 
and  possible  projects  for  the  future  with  several  groups  in  the  U.S.A.  on  such  and  related  research. 
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Electronic  absorption  spectra  of  B3  and  B3  in  neon  matrices 
and  ab  initio  analysis  of  the  vibronic  structure 

Muriel  Wyss,  Evgueni  Riaplov,  Anton  Batalov,  and  John  P.  Maiera) 
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Laboratoire  de  Chimie  Theorique,  Universite  de  Marne  la  Vallee,  F-77454  Champs  sur  Marne,  France 
(Received  17  July  2003;  accepted  4  August  2003) 

Mass  selected  /J3  ions  have  been  isolated  in  6  K  neon  matrices  and  their  absorption  spectra 
measured.  A  band  system  with  origin  at  467  nm  is  assigned  as  the  X  1A  [  electronic  transition 
of  the  cyclic  anion.  After  photobleaching,  the  1  2£'<—  X2A[  and  2  2E'  <—X2A[  band  systems  of 
neutral  cyclic  B3  are  observed  which  start  around  736  and  458  nm,  respectively.  Large  scale  ab 
initio  calculations  have  provided  potential  energy  surfaces  for  a  variational  treatment  of  the 
vibrational  motion.  Calculated  band  origins  leave  no  doubt  about  the  electronic  symmetry 
assignments.  The  complex  vibrational  structure  in  the  1  2E'  state,  which  is  due  to  relatively  strong 
Jahn-Teller  distortions,  appears  to  be  closely  reproduced  by  the  calculated  vibrational  energies  and 
intensities,  if  the  first  observed  stronger  line  is  identified  with  the  first  vibrationally  excited  state, 
placing  the  “true”  band  origin  of  the  1  2E'  state  at  775  nm  where  no  signal  with  significant  strength 
is  apparent.  The  2  2E'  state  undergoes  only  a  relatively  weak  Jahn-Teller  distortion  and  shows  a 
short  progression  with  an  observed  frequency  of  981(10)  cm-  that  compares  favorably  with  the 
theoretical  frequency  of  973  cm-1.  The  lE'  system  of  B3  shows  a  Jahn-Teller  activity  comparable 
to  that  of  the  1  2E'  state  of  B3 .  ©  2003  American  Institute  of  Physics.  [DOI:  10.1063/1.1613251] 


I.  INTRODUCTION 

There  has  been  a  growing  interest  in  the  structure  and 
energetics  of  small  boron  and  boron  containing  molecules 
because  they  have  a  variety  of  applications  such  as  for  the 
construction  of  high  temperature  semiconductor  devices,  as 
chemical  insulators,  and  in  the  production  of  refractory  ma¬ 
terials  and  explosives.1  However,  there  is  little  information 
on  the  spectroscopy  of  pure  boron  species  in  the  literature. 
Only  the  boron  dimer2  and  its  anion1  appear  to  be  well  char¬ 
acterized  experimentally  and  theoretically. 

Two  experimental  studies  have  been  carried  out  for  fi3, 
both  in  rare-gas  matrices.  One  of  these  is  an  electron  spin 
resonance  (ESR)  investigation  that  has  confirmed  the  pres¬ 
ence  of  three  equivalent  nuclei  in  the  2 A  [  ground  state4  and 
the  other  is  an  infrared  spectroscopic  study  in  which  the  de¬ 
generate  mode  (e')  has  been  observed  for  both  1 1 B 3  and 
(I/J3  isotopomers  (886  and  928  cm-1,  respectively).5  Several 
theoretical  investigations  on  the  fi3  molecule  have  predicted 
an  equilateral  cyclic  structure  (D3h)  with  a  2 ,4  J  electronic 
symmetry,6-9  and  the  vibrational  frequencies  in  the  elec¬ 
tronic  ground  state.  Of  the  low-lying  excited  electronic 
states,  only  the  states  of  2A3  and  2E'  symmetry  are  acces¬ 
sible  from  the  2 A  j  ground  state  by  dipole  allowed  excitation. 
Several  transition  bands  are  expected  in  the  near  infrared  and 
visible  regions.7,8  Because  the  Jahn-Teller  distortion  in  the  E 
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states  implies  stretching  of  a  strong  bond,  the  fi3  molecule 
should  provide  a  very  interesting  example  for  vibrational 
motion  that  proceeds  on  two  coupled  surfaces  even  for  the 
lowest  vibrational  quanta. 

Data  on  the  Bf  anion  are  even  more  scarce.  There  is 
only  a  theoretical  study  predicting  a  linear  'S.T  electronic 

°  Q 

ground  state  for  this  species  from  MP4  calculations.  How¬ 
ever,  MR-CI  calculations  carried  out  as  part  of  this  work 
show  clearly  that  the  lowest  energy  isomer  possesses  a  cyclic 
D3h  structure  with  a  *A]  electronic  ground  state.  As  far  as 
larger  negatively  charged  boron  molecules  are  concerned, 
only  the  gas  phase  photoelectron  spectrum  of  Bf  has  been 
reported.10 

In  the  present  work  the  electronic  absorption  spectra  of 
Bf  and  //3  have  been  measured  in  neon  matrices  at  6  K 
following  mass-selection.  In  order  to  provide  symmetry  as¬ 
signments  and  an  interpretation  of  the  complex  vibronic 
structure  observed  in  the  excited  electronic  state  1  2£’'  of  B  3 
as  result  of  the  Jahn-Teller  effect,  multireference  configura¬ 
tion  interaction  (MR-CI)  ab  initio  electronic  structure  calcu¬ 
lations  have  been  undertaken.  Based  on  diabatic  potential 
energy  surfaces  in  analytical  form,  the  vibrational  motion  has 
been  analyzed  in  variational  calculations  using  hyperspheri- 
cal  coordinates.  For  the  lower  transition  of  B3,  line  intensi¬ 
ties  and  a  tentative  vibronic  assignment  are  presented. 

II.  EXPERIMENT 

The  apparatus  employed  has  been  described 
previously.11  The  B  ,  anions  were  produced  in  a  cesium  sput- 

©  2003  American  Institute  of  Physics 
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FIG.  1.  Electronic  absorption  spectra  recorded  in  a  6  K  matrix  after  4  hours  of  mass-selected  codeposition  of  B}  with  neon.  The  bottom  trace  shows  the 
lE'— X  'A[  electronic  transition  of  ZJ3  overlapped  by  the  2  2E'—X2A J  system  of  B{ ,  produced  from  partial  neutralization  of  the  anions  impinging  on  the 
matrix  during  deposition.  The  top  trace  reveals  the  2  2E'-X2A  J  electronic  transition  of  B{  measured  after  exposure  to  UV  radiation.  Absorption  belonging 
to  the  anion  disappears. 


ter  ion  source  that  had  been  built  originally  for  the  produc¬ 
tion  of  carbon  anions.12  A  boron-carbide  (B4C)  rod  or 
pressed  pure  boron  powder  were  used  as  sputter  targets.  The 
ions  produced  in  the  source  were  guided  by  an  electrostatic 
lens  system  into  a  quadrupole  mass  analyzer,  where  the  B3 
anions  were  selected  and  then  codeposited  with  excess  of 
neon  onto  a  rhodium  coated  sapphire  plate  at  a  temperature 
of  6  K.  Ion  currents  of  25  nA  using  the  boron  sputter  target 
and  10  nA  from  B4C  could  be  obtained  for  the  mass  selected 
nB3  anions  with  unity  mass  resolution.  All  the  experiments 
but  one  were  performed  with  the  B4C  precursor  because  the 
available  pure  boron  sample  was  brittle. 

Matrices  of  about  150  /rm  thickness  were  grown  during 
4  hours.  Absorption  measurements  in  the  220-1100  nm 
range  were  carried  out  by  passing  monochromatic  light  from 
a  halogen  or  xenon  arc  lamp  through  the  matrix  employing  a 
waveguide-like  technique  and  achieving  a  path  length  of 
about  2  cm.  A  medium  pressure  mercury  lamp  (£^5.4  eV) 
was  used  to  neutralize  the  trapped  ions.  After  exposure  to  the 
UV  radiation  the  spectra  were  recorded  again. 

III.  OBSERVATIONS 
A.  Absorption  system 

Three  unknown  systems  built  upon  the  apparent  band 
origins  at  458,  467,  and  736  nm  were  observed  after  mass 


selected  (33  u)  codeposition  of  B3  ions  with  neon.  As  is 
illustrated  by  Fig.  1,  exposure  of  the  matrix  to  5.4  eV  radia¬ 
tion  resulted  in  the  complete  disappearance  of  the  467  nm 
system  and  the  growth  of  the  458  nm  one.  In  addition,  a  new 
system  starting  at  736  nm  evolved  (Fig.  2),  exhibiting  a  rich 
vibrational  structure  and  correlating  in  intensity  with  the  458 
nm  one.  Thus  one  concludes  that  the  carrier  of  the  467  nm 
system  is  the  B3  anion,  whereas  the  458  and  736  nm  band 
groups  are  due  to  the  neutral  B3  species.  Because  natural 
boron  contains  about  20  percent  as  its  10B  isotope,  it  was 
reasonable  to  check  whether  or  not  some  of  these  newly 
observed  absorptions  belong  to  the  B2C  molecule  or  its  ion. 
This  molecule  is  also  produced  in  the  cesium  sputter  source 
from  a  B4C  rod  and  its  10BnBC  isotopomer  has  the  same 
mass  as  B3  and  is,  therefore,  deposited  as  well.  Nevertheless, 
even  after  deposition  of  1 1 B2C  (34  u),  the  production  of 
which  is  much  higher  than  that  of  10BUBC-,  neither  of  the 
above  described  systems  were  seen  while  usage  of  a  pure 
boron  rod  in  the  source  led  again  to  the  detection  of  the  three 
band  systems  of  Figs.  1  and  2  with  the  same  intensity  ratios. 
This  supports  the  assignment  of  the  absorptions  observed  to 
the  B3  molecule  and  its  anion. 

The  equilateral  triangle  structure  of  the  neutral  molecule 
B3  in  the  ground  electronic  state  of  2 ,4  J  symmetry  is  well 
established.9  Dipole  transition  moments  of  symmetry  4" 
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FIG.  2.  Electronic  absorption  spectrum  of  the  1  2E'  —X  2A[  electronic  transition  of  B ,  recorded  after  4  hours  of  mass-selected  codeposition  with  neon 
followed  by  UV  irradiation  of  the  6  K  matrix. 


+E'  connect  this  state  only  with  states  of  symmetry  2 ,4  "  and 
2E' ,  respectively.  As  argued  in  detail  later  on,  only  2E' 
states  are  available  in  the  range  of  excitation  energies  inves¬ 
tigated  here. 

Because  the  cyclic  B3  neutrals  are  produced  by  photode¬ 
tachment  of  the  electron  from  the  mass-selected  B3  and 
isomerization  of  a  linear  to  cyclic  form  of  species  held  rela¬ 
tively  rigidly  in  a  6  K  neon  lattice  on  irradiation  with  3-5  eV 
photons  is  unlikely,  the  measured  electronic  spectmm  of  B3 
is  also  ascribed  to  a  cyclic  isomer  with  ground  state  of  ] 
symmetry.  This  is  in  line  with  the  observation  that  the  exci¬ 
tation  energy  is  close  to  the  corresponding  transition  in  neu¬ 
tral  B3,  as  to  be  expected  from  simple  molecular  orbital 
(MO)  considerations  since  the  same  MOs  are  involved  and 
the  theoretical  prediction  for  the  transition  (see  Sec.  IV). 
Therefore,  the  upper  state  of  the  band  system  observed  is 
assigned  to  the  1E'  electronic  symmetry. 

B.  Vibrational  structure 

The  vibrational  structure  apparent  in  the  observed  elec¬ 
tronic  spectra  (Figs.  1  and  2  and  Tables  I— III)  is  associated 
with  excitation  of  the  two  normal  modes,  ^(aj),  v2(e'),  in 
the  upper  electronic  states  because  at  6  K  the  population  is 
essentially  all  in  the  vibrationless  level  of  the  ground  state. 
Because  the  upper  electronic  states  of  fi3  and  fi3  have  E' 
symmetry  they  are  subject  to  Jahn-Teller  distortions  which 
may  cause  complex  vibrational  pattern. 

The  second  system  of  B3,  assigned  to  the 
2  2E' <—X2A[  transition,  has  the  simplest  pattern:  Starting 
with  the  origin  band  at  458  nm  a  seemingly  regular  progres¬ 
sion  mode  is  discernible  (Fig.  1)  with  two  spacings  that  are 


the  same  within  the  error  limit  (Table  III).  As  the  origin  band 
is  dominant,  no  large  geometry  change  is  expected  on  pass¬ 
ing  from  the  ground  to  the  excited  electronic  state. 

In  contrast,  the  vibronic  structure  of  the  first  system  of 
B3  (Fig.  1,  top  trace.  Table  II),  assigned  to  the  1  2E' 
<—  X  2 A  \  transition,  indicates  a  strong  Jahn-Teller  distortion: 
The  excitation  of  the  degenerate  v2(e')  Jahn-Teller  active 
mode  is  evident.  The  two  transition  systems  are  not  only 
very  different  in  complexity,  they  also  show  quite  different 
intensities:  The  spectra  indicate  an  absorbance  ratio  of  one 
order  of  magnitude  in  favor  of  the  second  band. 

The  absorption  spectmm  of  B3  (Fig.  2,  bottom  trace. 
Table  I),  assigned  to  the  1  1£,'<—  XlA[  transition,  shows  a 
complexity  similar  to  the  lower  transition  of  fi3.  The  pat¬ 
terns  of  these  two  systems  are  really  not  simple  and  a  sen¬ 
sible  assignment  can  only  be  made  by  means  of  a  theoretical 
treatment  as  attempted  in  Sec.  V. 


TABLE  I.  Positions  of  band  maxima  (X±0.2nm)  observed  in  the 
1  1Er -X  lA[  absorption  spectrum  of  B J  in  a  6  K  neon  matrix 


X.  (nm) 

v  (cm  *) 

Av  (cm  ') 

466.9 

21  412 

0 

448.4 

22  295 

883 

445.5 

22  440 

1028 

440.2 

22711 

1299 

436.2 

22  919 

1507 

431.9 

23  147 

1735 

428.4 

23  336 

1924 

427.1 

23  407 

1995 

418.4 

23  894 

2482 
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TABLE  II.  Positions  of  band  maxima  (\±0.2nm)  in  the  1  2E'  —X2A[ 
absorption  spectrum  of  S3 ,  vibronic  symmetry  I  and  calculated  intensities. 


\(nm) 

v  (cm  *) 

AiXcm  1 

Experiment11 

) 

Theory 

r 

Intensity 

0 

E' 

1.00 

230 

K 

735.8 

13  587 

0 

685 

685 

E' 

1.45 

733 

A' 

715.5 

13  972 

385 

1070 

1074 

E' 

0.74 

1073 

A[ 

705.6 

14  168 

581 

1266 

1269 

E' 

0.59 

700.7 

14  268 

681 

1366 

1337 

E' 

1.80 

1425 

a; 

1660 

A'i 

681.6 

14  667 

1080 

1765 

1732 

E' 

0.25 

678.3 

14  739 

1152 

1837 

1818 

E' 

1.70 

1801 

K 

1910 

A'i 

673.4 

14  846 

1259 

1944 

1933 

E' 

0.60 

2123 

E' 

0.20 

2195 

A[ 

2202 

E' 

0.01 

2310 

E' 

0.00 

2340 

a'2 

2384 

a; 

650.0 

15  380 

1793 

2478 

2444 

E' 

0.98 

2481 

E' 

0.01 

2537 

a; 

2564 

a; 

2648 

A) 

2714 

A2 

2728 

E' 

0.03 

2762 

a; 

2834 

E' 

0.11 

2848 

E' 

0.04 

2868 

^2 

2915 

a; 

2920 

E' 

0.45 

3057 

E' 

0.02 

“For  two  different  band  origins,  see  text. 


IV.  THEORY 

Symmetry  assignments  and  analysis  of  complex  vibra¬ 
tional  structures  as  observed  for  B3  require  the  knowledge  of 
sufficiently  accurate  and  extended  potential  energy  surfaces 
(PES)  for  low-lying  electronic  states.  Because  previous  the¬ 
oretical  efforts  have  been  limited  to  a  few  cuts  through  se¬ 
lected  PESs,7  we  have  performed  large-scale  ab  initio  calcu¬ 
lations  for  the  states  relevant  to  the  absorption  spectra 
presented  here.13 

The  orbital  ordering  for  B3  in  D3h  symmetry  is  laj, 
le' ,  2a[,  2e',  la",  3 a\  for  the  orbitals  occupied  in  the 
ground  state,  and  3e',  le",  4er,  1  a'2,...  for  the  lower  vir¬ 
tual  orbitals.  Clearly,  the  ground  state  is  of  2  A  [  symmetry 
and  excited  states  accessible  by  dipole  selection  rules  have 
2 A"  and  2E'  symmetries.  While  the  1  2 A"  state  is  dominated 
by  the  single  configuration  resulting  from  the  3a  la2  ex¬ 
citation,  two  low-lying  2E'  states  should  result  from  substan¬ 
tial  mixing  of  the  configurations  corresponding  to  the  exci¬ 
tations  3e'<— 3a]  and  3a]<— 2e',  respectively.  The  ground 
state  configuration  of  B3  is  that  of  B3  but  with  an  extra 
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TABLE  III.  Positions  of  band  maxima  (A.  ±  0.2  nm )  observed  in  the 
2  2E'-X2A[  absorption  spectrum  of  B 3 .  calculated  vibrational  energies 
and  vibronic  symmetry  L. 


A.(nm) 

v  (cm  J) 

Aj'fcm  *) 

Expt. 

Theory 

r 

458.0 

21828 

0 

0 

E' 

848 

^2 

438.3 

22809 

981 

973 

E' 

974 

a; 

1258 

E' 

1746 

E' 

1815 

^2 

420.1 

23797 

1969 

1951 

E' 

2018 

a; 

2144 

E' 

electron  in  the  highest  occupied  orbital  3a  j .  The  lowest  di¬ 
pole  allowed  excitation  is  here  dominated  by  the  3e'<— 3a] 
promotion. 

A  configuration  mixing  as  indicated  above,  on  top  of  the 
interaction  between  the  two  components  of  the  2E '  states  for 
distorted  C s  geometries  where  both  have  2 A  ’  local  symme¬ 
try,  restricts  our  choice  of  method  effectively  to  the  varia¬ 
tional  multireference  configuration  interaction  (MR-CI)  pro¬ 
cedure.  For  efficiency  we  use  a  complete  active  space  (CAS- 
SCF)  for  the  reference  wave  functions  and  the  internally 
contracted  singles  and  doubles  Cl,  as  implemented  in  the 
MOLPRO14  program  package.  The  whole  set  of  geometries  is 
treated  with  an  active  space  spanned  by  the  orbital  set  from 
2a  J  to  le”;  extended  calculations  for  a  subset  of  C2v  geom¬ 
etries  use  an  active  space  that  includes  the  4e'  and  1  a'2 
orbitals  as  well.  The  basis  set  is  constructed  from  the  aug¬ 
mented  correlation-consistent  valence  triple  zeta  (aug-cc- 
pVTZ)  set  for  the  B  atom,15  further  augmented  by  two  s 
functions  (exponents  0.7  and  0.23)  and  one  p  set  (exponent 
0.35)  at  the  bond  centers.  Comparison  between  experiment 
and  theory  for  B2  suggests  that  this  computational  procedure 
supports  an  accuracy  of  about  0.02  bohr  in  bond  lengths  and 
about  20  cm-1  in  vibrational  frequencies.  These  calculations 
have  been  performed  for  a  grid  of  geometries  which  was 
defined  by  bond  length  variations  in  steps  of  0.05  bohr,  start¬ 
ing  from  the  reference  bond  length  of  r  =  3.00  bohr  and  ex¬ 
tending  to  all  geometries  for  which  one  of  the  relevant 
electronic  energies  was  less  than  3000  cm-1  above  the  cor¬ 
responding  reference  energy. 

The  adiabatic  PES’s  have  been  cast  in  analytical  form  by 
least-squares  fits  to  polynomials  in  the  symmetric  stretch  and 
radial  deformation  coordinates  and  harmonic  functions  of  the 
azimuthal  deformation  coordinate,  as  described  in  previous 
work  on  alkali  trimers.16'17  In  contrast  to  the  alkali  trimers, 
Jahn-Teller  distortion  in  the  lower  E'  states  of  B3  and  B3  is 
weak  as  compared  to  vibrational  excitation,  so  that  both 
sheets  of  the  PES  for  degenerate  E  states  are  substantially 
populated  even  for  low  vibrational  energies.  In  order  to  avoid 
the  notorious  complications  caused  by  the  geometric  phase 
in  an  adiabatic  treatment,  we  use  a  diabatic  framework  for 
the  vibrational  motion,  i.e.,  diabatic  surfaces  are  constructed 
from  the  adiabatic  ones  by  means  of  2X2  rotations  with 
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angles  cf>l2,  where  </)  is  the  pseudo-rotational  coordinate.  Hy- 
perspherical  coordinates  are  still  the  proper  choice  for  a  cy¬ 
clic  trimer  but  the  computational  procedures  used 
previously18  had  to  be  extended  to  account  for  pairs  of 
coupled  diabatic  surfaces.  Both,  the  analytical  fits  and  the 
vibrational  calculations  introduce  numerical  errors  in  the  or¬ 
der  of  only  1  cm-1  for  the  first  transition  in  B3,  i.e.,  consid¬ 
erably  smaller  errors  than  those  from  the  electronic  structure 
calculations.  For  the  other  upper  states  our  fits  with  standard 
errors  of  about  10  cm-1  are  only  preliminary  because  of 
complications  from  additional  conical  intersections  and/or 
notorious  intruder  state  problems  which  cause  unphysical 
humps  in  the  upper  state  surfaces  in  particular  for  B ,  . 

Vibrational  line  intensities  are  calculated  from  transition 
dipole  moment  surfaces  which  are  again  expanded  in  poly¬ 
nomials  as  outlined  above  for  the  PES.  Because  the  observed 
bands  stem  from  in-plane  transition  moments  (parallel 
bands),  vibronic  selection  rules  leave  only  2E'  vibronic 
states  accessible  in  the  case  of  free  B2 .  The  degree  of  viola¬ 
tion  of  this  selection  rule  due  to  matrix  effects  is  difficult  to 
assess. 

A  sensible  assignment  of  vibronic  states  by  some  kind  of 
vibrational  quantum  numbers  is  difficult  in  cases  of  strong 
vibronic  coupling,  i.e.,  if  linear  and  quadratic  Jahn-Teller 
(JT)  terms  are  of  comparable  size.  Reference  can  be  made  to 
two  simple  limiting  cases:  (i)  Strong  quadratic  terms  produce 
three  equivalent  minima  on  the  lower  (adiabatic)  surface 
which  are  separated  by  a  relatively  high  barrier.  The  vibra¬ 
tional  wave  functions  for  the  lower  sheet  can  then  be  consid¬ 
ered  a  linear  combination  of  three  equivalent  local  vibrators 
with  C2v  symmetry  placed  at  the  respective  minima  and  can 
be  labeled  by  C2v  quantum  numbers  u  , ,  v2 ,  and  v3  for 
symmetric  stretch,  symmetric  bend  and  asymmetric  stretch, 
respectively,  (ii)  Weak  quadratic  JT  terms  produce  only  neg¬ 
ligible  localizing  barriers  for  the  pseudorotational  motion. 
Then  the  quantum  numbers  for  a  free  internal  rotation  are 
adequate,  i.e.,  vt,  v2  for  symmetric  stretch  and  deformation 
and  1 2  for  the  pseudorotation  angular  momentum.  In  an  adia¬ 
batic  frame,  v2  and  l2  are  half  integer  numbers  due  to  the 
geometric  phase  in  electronic  and  vibrational  factors  of  the 
vibronic  wave  function.  The  finger  print  for  case  (i)  is  a  near 
degeneracy  of  E-A  pairs,  that  for  case  (ii)  is  a  near  degen¬ 
eracy  of  A  ]  —  A  2  pairs.  The  actual  assignment  of  model  quan¬ 
tum  numbers  is  made  by  inspection  of  plots  of  the  vibra¬ 
tional  wave  functions.  The  model  following  case  (i)  has  been 
successfully  used  for  rovibronic  analysis  of  alkali  trimer  E" 
<—E'  transition  bands;16-18  a  thorough  discussion  of  this 
model  has  been  given  there.  Further  details  of  our  computa¬ 
tional  procedures  may  be  found  elsewhere.13 

V.  ASSIGNMENT  AND  DISCUSSION 

Our  CAS-SCF  calculations  of  vertical  electronic  transi¬ 
tion  energies  at  the  i)3h  reference  geometry  (bond  length 
3.00  bohr)  reproduce  the  ordering  of  states  as  previously 
obtained,  also  from  CAS-SCF:9  X2A[,  1  2A2,  1  2E',  1  4E", 
1  4£\  1  2E",  2  2E' ,  1  4A[ ,  1  4A2  ...  .  Pertinent  data  of  the 
PESs,  as  obtained  in  the  MR-CI  calculations  for  the  states 
relevant  here,  are  collected  in  Table  IV.  Figure  3  shows  cuts 


TABLE  IV.  Transition  energies,  Jahn-Teller  distortion  energies  and  equilib¬ 
rium  geometries  for  selected  states  of  B3  from  MR-CI  with  the  large  active 
space  (energies  in  cm-1,  distances  in  bohr  and  angles  in  degree). 


State: 

i2a; 

1  2E' 

2  2E' 

Tv 

13  781 

22  291 

T„ 

13  057 

22  202 

Too 

12  829 

22  272 

ZPE 

1457 

1207 

1522 

^■stab 

475 

39 

£loc 

187 

9 

rmin 

2.968 

2.982 

3.015 

0mm 

60.0 

65.3 

61.9 

Gadd 

3.093 

2.982 

^sadd 

57.2 

60.9 

of  these  potentials  for  the  bending  deformation  mode,  which 
illustrate  in  particular  the  Jahn-Teller  distortions  in  the  E' 
states.  The  JT  stabilization  in  the  1  2E'  state  of  475  cm-1  is 
on  the  order  of  stabilizations  observed  for  alkali  trimer  2E" 
states,  but  it  amounts  to  little  more  than  1/3  of  the  zero  point 
energy  and  only  about  2/3  of  the  deformation  mode  energy 
quantum.  The  pseudorotational  barrier  between  the  three 
equivalent  minima  is  only  187  cm-1  but  it  is  sufficiently 
broad  to  effectively  lift  the  AXIA2  degeneracy,  indicating  a 
case  (i)  description  of  the  vibronic  states.  In  contrast,  for  the 
2  2E'  state  with  stabilization  and  localization  energies  of 
only  39  and  9  cm-1,  respectively,  B3  can  be  considered  as 

Angle 


Symmetric  bend  :  (2r3-  r2-  r,)  [a0] 

FIG.  3.  One-dimensional  cuts  of  the  PESs  along  the  bend  coordinate  for  C2v 
geometries  with  (r1  +  r2  +  r3)/3  =  3.00  bohr  for  the  X2A[  and  the  two  2E' 
electronic  states  of  B3  . 
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FIG.  4.  Comparison  of  the  calculated  positions  and  intensities  of  the  Jahn- 
Teller  coupled  vibronic  levels  for  the  1  2E'  state  of  U  t  with  the  matrix 
spectrum. 


nearly  freely  pseudo-rotating,  indicating  a  case  (ii)  descrip¬ 
tion.  This  state  has  a  second  conical  intersection  of  its  two 
electronic  components  at  C2v  geometries  with  apex  angles 
around  68°.  Such  a  feature  causes  severe  complications  for 
the  fitting  of  a  PES.  Our  present  analytical  potential  is  reli¬ 
able  only  below  that  crossing  and  furnishes  but  a  few  low- 
lying  vibrational  levels.  The  shapes  of  these  PES,  and  even 
more  so  the  rather  unexpected  shapes  of  the  transition  mo¬ 
ment  surfaces,  are  consequences  of  the  strong  mixing  of  the 
two  configurations  that  are  formed  by  the  3e'  <—3a[  and  the 
3a[<—  2e'  excitations.  Actually,  at  the  D 2h  reference  geom¬ 
etry  the  mixing  ratios  are  62:58  and  61:59  for  the  1  2E'  and 
2  2E '  states,  respectively.  The  corresponding  dipole  transi¬ 
tion  moments  of  0.137  and  0.253  a.u.,  respectively,  indicate  a 
clear  difference  in  the  intensities  of  the  two  bands,  as  is 
indeed  observed. 

The  calculated  positions  of  the  transition  bands  are  cru¬ 
cial  for  the  assignment  of  electronic  symmetries.  We  find  the 
1  2 E'  band  origin,  7 0(l ,  at  12  829  cm-1  from  the  MR-CI  with 
extended  reference  space.  The  corresponding  results  from  the 
smaller  MR-CI  and  from  a  Davidson  correction  are  12  988 
and  12  728  cm-1,  respectively,  indicating  an  uncertainty  in 
Too  of  about  250  cm  1 .  However,  the  observed  spectrum 
starts  with  a  line  at  13  587  cm-1,  which  is  752  cm-1  higher 
than  calculated.  This  difference  is  too  large  to  be  attributed  to 
matrix  effects.  It  is  not  only  larger  than  the  estimated  uncer¬ 
tainty,  but  its  sign  is  unlikely  since  in  variational  calculations 
excited-state  energies  are  usually  less  well  accounted  for 
than  ground-state  energies.  This  suggests  that  the  first  ob¬ 
served  line  is  not  the  band  origin  but  corresponds  to  the  first 
excited  vibrational  level.  This  is  indeed  supported  by  the 
positions  and  intensities  calculated  for  the  higher  vibrational 
states,  as  obvious  from  Fig.  4,  where  the  observed  spectrum 
is  compared  with  the  theoretical  one,  the  latter  shifted  only 
by  30  cm-1.  Agreement  appears  to  be  rather  satisfactory  con¬ 
sidering  the  well-known  fact  that  matrix  effects  may  affect 
the  band  positions  and  in  particular  their  intensities.  The  trip¬ 
let  feature  seen  in  the  center  of  the  spectrum  belongs  to  vi¬ 
brational  states  that  are  heavily  mixed,  so  that  the  calculated 
intensities  may  well  need  some  redistribution,  since  the  com¬ 


position  of  the  vibronic  wave  functions  is  computationally 
more  sensitive  to  the  accurate  description  of  the  coupling 
than  the  vibronic  energies.  However,  this  assignment  of  the 
vibrational  structure  is  in  conflict  with  the  fact  that  the  ob¬ 
served  spectrum  does  not  show  a  line  with  significant 
strength  at  the  predicted  position  of  the  band  origin.  On  the 
other  hand,  if  the  band  origin  is  identified  with  the  first 
strong  observed  line,  we  can  not  explain  that  two  strong 
observed  lines  do  not  have  theoretical  counterparts  and  that  a 
line  with  a  strong  calculated  intensity  is  missing.  It  is  diffi¬ 
cult  to  see  how  matrix  effects  could  be  held  responsible  for 
this  mismatch.  A  gas  spectrum  is  necessary. 

Table  II  lists  observed  and  calculated  vibrational  ener¬ 
gies  and  the  calculated  relative  intensities,  making  also  use 
of  the  theoretical  first  excitation  energy  to  define  the  “true” 
band  origin,  which  becomes  12  902  cm-1.  Figure  5  displays 
vibrational  levels  ordered  according  to  vibronic  symmetry 
and  assigned  with  vibrational  quantum  numbers  of  a  C2u 
vibrator  as  outlined  above.  This  assignment  is  reasonable  for 
the  lower  levels  but  loses  meaning  for  higher  ones,  in  par¬ 
ticular  if  levels  are  close  and  undergo  substantial  mixing. 
This  is  regularly  the  case  for  neighboring  levels  of  polyads 
formed  by  replacing  two  quanta  of  symmetric  bend  by  one 
quantum  of  asymmetric  stretch.  Note  that  only  the  levels  of 
vibronic  E  character  are  accessible. 

The  calculated  band  origin  for  the  2  2E '  state  of  22  272 
cm-1  agrees  fairly  well  with  the  observed  position  of  the 
main  line,  21  828  cm-1.  There  is  no  other  candidate  close  by. 
Vibrational  energy  levels,  which  are  only  preliminary  in  this 
case  due  to  problems  with  the  fit  of  the  PES,  are  included  in 
Table  III.  The  observed  line  positions  may  be  nicely  lined  up 
with  calculated  values  but  this  is  not  really  definitive  as  long 
as  theoretical  intensities  are  lacking  for  this  state.  While  the 
ground  vibrational  state  is  93%  a  [ ,  the  first  upper  state  is 
60%  e[  and  30%  a  [ ,  i.e.,  the  progression  observed  involves 
mainly  the  deformation  mode  with  only  zero-point  pseudoro¬ 
tation.  The  proper  assignment  according  to  case  (ii)  is  (0  1/2 
1/2),  (0  3/2  1/2),  and  (0  5/2  1/2)  for  the  three  levels  involved. 

For  the  1  2A2  state  we  calculate  a  vertical  excitation  en- 
ergy  Tv  of  6300  cm-1.  This  region  was  inaccessible  in  the 
experiments  carried  out  hitherto.  As  a  general  comment  on 
the  transition  energies,  we  note  that  dynamical  electron  cor¬ 
relation,  accounted  for  only  by  the  singles-doubles  Cl  build 
on  top  of  the  CAS-SCF,  causes  changes  of  up  to  2000  cm-1. 
Comparison  to  previously  published  theoretical  excitation 
energies,  which  are  limited  to  CAS-SCF,  shows  indeed  dif¬ 
ferences  of  this  order.8 

In  the  case  of  the  B3  anion  the  calculated  adiabatic  ex¬ 
citation  energy  Te  of  23  089  cm~ 1  differs  from  the  observed 
band  origin  at  21413  cm-1  by  1676  cm-1.  Less  than  300 
cm-  may  be  attributed  to  the  change  in  the  zero  point  en¬ 
ergy  (ZPE),  leaving  a  defect  of  around  1400  cm-1.  Extension 
of  our  basis  set  by  diffuse  functions,  which  might  be  impor¬ 
tant  for  an  excited  state  of  an  anion,  reduces  this  defect  by 
only  about  100  cm-1  but  the  Davidson  correction  yields  a 
surprisingly  large  further  reduction  of  about  800  cm~ 1  in  this 
case.  The  assignment  to  the  1  X  lA[  electronic  transi¬ 

tion  of  cyclic  B 3  is  not  questioned  since  no  other  candidates 
are  available.  A  vibrational  analysis  of  this  transition  has  to 
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FIG.  5.  Pattern  of  Jahn-Teller  coupled  vibronic  levels  of  the  1  2E'  state  in 
B3 .  For  the  definition  of  the  modes  see  text. 

be  deferred  until  a  better  analytic  representation  of  the  upper 
state  PES  is  available. 

VI.  CONCLUSION 

The  matrix  isolation  technique  in  combination  with  mass 
selection  has  been  applied  to  observe  for  the  first  time  the 
electronic  spectra  of  the  boron  trimer  and  its  anion.  As  a 
result,  the  l2E'<—  X2A[  and  the  2:E'<—  X2A[  electronic 
transitions  for  B 3  and  *£'<—  X  lA[  for  Z?3  have  been  identi¬ 
fied.  The  vibronic  structure  of  the  1  2E'  electronic  state  of 
the  neutral  species  and  the  1  lE'  state  of  the  anion  exhibit  a 
strong  influence  of  the  Jahn-Teller  effect.  Ab  initio  calcula¬ 
tions  of  potential  energy  surfaces  for  four  electronic  states  of 
fi3  and  two  electronic  states  of  B,  have  validated  the  assign¬ 
ment  of  the  observed  bands  and  provided  the  basis  for  an 


analysis  of  the  vibrational  structure  of  the  transition  to  the 

I  2E'  state  of  a  free  fi3.  The  observed  pattern  can  be  con¬ 
sidered  as  well  reproduced  if  the  true  origin  of  the  transition 
is  placed  at  about  12  900  cm-1  where  no  significant  absorp¬ 
tion  line  is  seen.  An  in-depth  analysis  will  require  high  reso¬ 
lution  gas-phase  measurements  and  theory  developments. 
The  present  work  provides  the  necessary  basis  for  these  to  be 
undertaken. 
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The  rotationally  resolved  origin  band  in  the  2  2E’  <—  X  2A  [  electronic  spectrum  of  cyclic  B3  has  been 
observed  by  cavity  ring  down  spectroscopy  in  the  gas  phase.  The  B3  molecule  was  generated  in  a 
supersonic  planar  plasma  containing  decaborane  (B10H14)  and  neon  as  a  carrier  gas.  The  rotational 
structure  pattern  is  that  of  a  cyclic  molecule.  It  is  analyzed  assuming  an  equilateral  triangle  in  both 
electronic  states.  The  band  origin  is  determined  to  be  21  853.52  cm~  ,  and  the  bond  lengths 
1.603  77(106)  A  in  the  ground  and  1.619  07(96)  A  in  the  excited  electronic  state  are  inferred  from 
analysis  of  the  rotational  structure.  ©  2004  American  Institute  of  Physics. 
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I.  INTRODUCTION 

Well  resolved  absorption  and  emission  spectra  of  cyclic 
trimers  have  been  reported  only  for  few  species,  including, 
Al3 , 1  Li3,2  Cu3,3  Hi  ,4  and  Na3.5  As  the  small  boron  mol¬ 
ecules  have  applications  in  industry  as  chemical  insulators, 
explosives,  semiconductors,  and  high  energy  density  materi¬ 
als,  their  structure  is  of  relevance.6'7  The  boron  trimer  has 
been  investigated  only  in  rare  gas  matrices.  First  came  the 
electron-spin-resonance  observation  which  showed  that  B3 
contains  an  unpaired  electron  and  that  the  three  nuclei  are 
equivalent.8  The  infrared  spectrum  of  B3  was  also  observed,1 
and  recently  the  2  2E' *—X2A[  electronic  transition  of  B3  in 
a  neon  matrix  was  identified  using  a  mass  selective 
approach.9  However,  there  is  no  spectroscopic  information 
on  the  boron  trimer  B3  in  the  gas  phase. 

There  is  interest  in  the  spectrum  of  cyclic  B3  because 
several  of  the  excited  electronic  states  are  degenerate  in  D3h 
symmetry.  According  to  the  Jahn-Teller  theorem,  there  exists 
at  least  one  vibrational  mode  that  removes  the  electronic  de¬ 
generacy  and  results  in  the  stabilization  of  the  lowest  com¬ 
ponent  by  lowering  the  symmetry  of  the  equilibrium  struc¬ 
ture.  Most  recent  theoretical  calculations  on  B3  predicted  an 
equilateral  structure  (D3h)  with  an  2A[  electronic  symmetry 
in  the  electronic  ground  state  and  an  obtuse  triangle  (C2u) 
with  an  apex  angle  3*60°  for  the  2  2E '  excited  electronic 
state.810'11  However,  a  low  stabilization  energy  (difference 
between  the  conical  intersection  and  the  minima)  compared 
with  zero  point  energy  for  the  2  2E'  state  is  reported.9  There¬ 
fore  the  B3  molecule  can  be  considered  to  have  D3h 
symmetry  in  the  2  2E'  state  as  an  effective  structure  and  the 
2  2E' <—  X 2 A  [  electronic  system  can  be  considered  as  a  quasi 
“ D3h-D3h ”  transition. 

In  the  present  work  we  report  the  observation  of  the 
rotationally  resolved  origin  band  of  the  2  2£',<—  X2A[  elec¬ 
tronic  spectrum  of  B3 .  The  analysis  of  the  spectrum  gives  an 
estimate  of  the  molecular  geometry;  the  first  gas  phase  struc¬ 
tural  determination. 


“'Author  to  whom  correspondence  should  be  addressed;  Fax:  +41-61- 
2673855.  Electronic  mail:  J.P.Maier@unibas.ch 


II.  EXPERIMENT 

This  consists  of  a  standard  cavity  ring  down  setup  sam¬ 
pling  a  plasma  generated  in  a  pulsed  supersonic  slit  jet 
expansion.12  A  small  fraction  of  pulsed  laser  light  from  a 
tunable  dye  laser  is  coupled  into  an  optical  cavity  consisting 
of  two  highly  reflective  mirrors  (R> 99.99%).  The  fraction 
of  light  leaking  out  of  the  cavity  has  an  envelope  which  is  a 
first-order  exponential  decay,  exp(— t/f).  The  ring  down  time 
t  reflects  the  absorption  coefficient.  A  spectrum  is  recorded 
by  measuring  t  as  function  of  the  wavelength  of  the  laser 
(running  at  30  Hz).  Forty-five  ring  down  events  are  averaged 
at  each  wavelength  before  the  digitized  data  are  downloaded 
to  a  workstation.  The  —0.05  cm-1  linewidth  of  the  laser  was 
attained  with  an  internal  etalon. 

The  pulsed  slit  jet  system  is  located  in  a  vacuum  cham¬ 
ber  evacuated  by  a  roots  blower.  The  plasma  is  produced  in  a 
nozzle  incorporating  a  discharge  in  the  high  pressure  expan¬ 
sion.  The  orifice  comprises  a  metal  plate  anode  and  two 
sharp  stainless  steel  cathodes  that  form  the  actual  slit  (30 
mmX300  /rm).  A  pulsed  negative  voltage  (—800  V/100  mA 
and  100  /rs  in  duration)  is  applied.  The  geometry  of  the 
orifice  is  such  that  the  discharge  is  confined  upstream  of  the 
supersonic  expansion.  To  produce  B3,  solid  B10H14  is  heated 
in  the  oven  to  70±10°C  with  neon  as  a  carrier  gas  at  a 
backing  pressure  of  10  bars.  In  order  to  avoid  condensation 
inside  the  source,  the  valve  and  body  of  the  nozzle  are  heated 
as  shown  in  Fig.  1.  The  B3  molecules  are  efficiently  pro¬ 
duced  at  a  distance  5  mm  down  stream  from  the  slit.  The 
spectrum  is  calibrated  via  the  neon  atomic  lines  that  are  ob¬ 
served  in  the  jet  expansion. 

III.  RESULTS 

The  2  2E'*—  X  2 A  [  electronic  spectrum  of  B3  was  ini¬ 
tially  identified  in  absorption  in  a  6  K  neon  matrix  using 
mass  selection  with  origin  band  located  at  21  828  cm-19 
Based  on  this  the  band  was  sought  in  the  gas  phase.  The 
absorption  signal  was  subsequently  found  near  21  900  cm-1 
(Fig.  2)  and  the  observed  profile  of  the  band,  —25  cm-  in 
width,  was  in  a  good  agreement  with  preliminary  simula¬ 
tions.  No  other  signal  could  be  detected  within  ±300  cm-1 
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Plasma 

FIG.  1.  Schematic  diagram  of  the  source  used  to  produce  B3.  Decaborane 
was  heated  up  to  70±10  °C  in  the  oven.  Neon  at  a  backing  pressure  of  10 


bars  was  used  as  a  carrier  gas.  The  pulsed  plasma  was  produced  by  a  slit 
nozzle  to  which  the  high  voltage  was  applied. 


from  the  observed  band.  Though  the  spectrum  visually  ap¬ 
pears  “noisy”  (typical  S/N  of  at  least  5)  all  the  peaks  are 
genuine.  This  is  shown  by  the  top  trace  which  is  an  indepen¬ 
dent  recording  of  this  part  of  the  band.  Corresponding  traces 
were  obtained  for  other  regions  proving  that  all  the  lines 
given  in  Tables  I  and  III  are  real. 

The  spectrum  exhibits  rotationally  resolved  P  and  R 
branches  slightly  blended  with  strong  neon  atomic  lines.  The 
Q  branch  is  not  visible  due  to  a  complex  splitting  and  blend¬ 
ing  with  an  intense  neon  line.  Actually,  it  is  a  composition  of 
four  spectra,  as  two  boron  isotopes  are  present.  Because  the 
natural  isotopic  abundance  of  nB:  10B  in  the  precursor  is  4:1, 
the  relative  concentration  of  the  two  isotopomers  33B3 : 32B3 
(two  equivalent  structures  of  10BnB2)  should  be  4:3.  The 
peaks  in  the  spectra  of  31B3  and  30B3  could  not  be  assigned 
due  to  their  small  abundance.  As  a  result,  the  observed  spec¬ 
trum  is  essentially  an  overlap  of  the  two  isotopomers  33B3 
and  32B3 . 

The  B3  cyclic  molecule  in  the  excited  degenerate  state 
should  undergo  a  deformation  due  to  the  Jahn-Teller  effect, 
however,  the  distortion  from  an  equilateral  triangle  can  be 
neglected  as  a  result  of  low  stabilization  energy  [merely  39 
cm-  compared  with  a  zero-point  energy  (1522  cm-1)  in  the 
2  2E"  state9].  A  similar  situation  was  encountered  for  the 
cyclic  C3H  radical,  where  deformations  of  the  C2v  structure 
were  neglected  due  to  small  energy  differences  between  the 
Cs  and  C2v  forms.13  Therefore,  it  is  reasonable  to  assume  that 
the  effective  symmetry  of  the  2  2E'  state  is  D 3 h  for  the 
analysis  of  the  rotational  structure  in  the  observed  spectrum. 


FIG.  2.  Rotationally  resolved  2  2E'  —  X  2A[  electronic  origin  band  of  B3  in 
the  gas  phase  measured  by  cavity  ring  down  spectroscopy  through  a  super¬ 
sonic  planar  plasma.  The  two  upper  traces  are  independent  recordings  to 
show  that  the  peaks  are  real  and  not  noise.  The  three  bottom  traces  are  the 
simulated  patterns  33B3  +  32B3  (in  ratio  4:3).  32B3  and  33B3 .  The  two  strong 
lines  off  scale  belong  to  neon.  Bands  marked  with*  are  unassigned.  Lines  in 
the  simulation  for  which  counterparts  in  the  spectrum  could  not  be  found  are 
marked  with  “  + 

The  measured  positions  of  the  33B3  lines  are  given  in 
Table  I.  The  spectral  line  frequencies  were  analyzed  using  a 
conventional  Hamiltonian  for  a  symmetric  top  molecule.  In 
order  to  analyze  the  imperfectly  resolved  rotational  spectrum 
it  was  assumed  that  the  structure  of  the  B3  molecule  has  I)2h 
symmetry  in  both  states.  This  way  the  number  of  spectro- 


TABLE  I.  Frequencies  (in  cm  ')  and  assignments  in  the  2  2£“<—  X  2 A  [ 
absorption  spectrum  of  33B3 . 


nk 

^obs 

^calc 

^obs  ^calc 

1  r22 

21  850.374 

21  850.543 

-0.169“ 

22-33 

21  849.223 

21  849.296 

-0.073 

33-44 

21  848.070 

21  848.027 

0.043 

44-55 

21  846.731 

21  846.736 

-0.005 

66-7  7 

21  844.118 

21  844.086 

0.032 

77-88 

21  842.710 

21  842.727 

-0.017 

IfOo 

21  855.230 

21  855.305 

-0.075 

22-l  i 

21  856.499 

21  856.440 

0.059 

33-22 

21  857.600 

21  857.552 

0.048 

44-33 

21  858.629 

21  858.642 

-0.013 

55-44 

21  859.690 

21  859.710 

-0.020 

88-77 

21  862.770 

21  862.778 

-0.008 

kO 

00 

21  863.713 

21  863.756 

-0.043 

1010-99 

21  864.749 

21  864.712 

0.037 

“Not  included  in  the  least-squares  fit. 
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TABLE  II.  Bond  lengths  and  spectroscopic  constants  of  the  33B3  molecule 
inferred  from  the  partial  rotational  analysis  of  the  2  2E '  <—  X  2  A  [  electronic 
origin  band.” 


2a; 

2  2E' 

Obs. 

Calc.b 

Obs. 

Calc.b 

r 0  (A) 

1.603  77(106) 

1.619  07(96) 

d0  (deg) 
re  (A) 

60. 0C 

1.5706 

60.0“ 

1.5955 

de  (deg) 

60.0 

60.8 

Too  (cm-1) 

21  853.52d 

22  723 

B  (cm-1) 

1.190  64(157)“ 

1.168  25(138)“ 

C  (cm-1) 

0.595  32(79)“ 

0.584  12(69)“ 

“The  numbers  in  parentheses  represent  the  one  standard  deviation  derived 
from  the  least  square  fitting.  A  rms  of  the  fitting  is  0.043  cm  . 
bReference  9. 

“Fixed. 

dError  of  calibration  is  around  0.01  cm-1. 

“The  rotational  constants  were  calculated  from  the  determined  molecular 
structures. 


scopic  parameters  could  be  reduced  to  3;  the  bond  lengths 
r',  r"  and  band  origin  7 0(l  were  then  determined  directly 
from  a  least-squares  fit  to  the  line  positions. 

The  results  for  the  33B3  molecule  are  listed  in  Table  II. 
The  determined  bond  lengths  (r0)  are  2%  longer  than  the 
calculated  ones  (rc)  (Ref.  9)  in  both  the  ground  and  excited 
states.  In  the  ground  state  this  is  because  of  zero-point  vibra¬ 
tions  as  the  analysis  assumes  a  negligible  inertial  defect.  The 
difference  in  the  excited  state  is  additionally  caused  by  im¬ 
perfections  of  the  simple  model  used  (assuming  an  equilat¬ 
eral  triangle  0=60°  and  ignoring  the  lahn-Teller  effect)  and 
anharmonicity  of  the  potential. 

The  observed  rotational  profile  was  also  simulated  with 
the  obtained  molecular  constants  using  the  program  WANG. 14 
The  simulated  spectrum  of  l3B3  is  shown  in  Fig.  2.  Best 
agreement  with  the  experiment  was  attained  with  a  rotational 
temperature  of  30  K.  A  linewidth  of  0.05  cm-1  was  used  for 
the  simulation,  though  in  the  observed  spectrum  it  is  larger 
due  to  Doppler  and/or  lifetime  broadening. 

In  the  case  of  the  32B3  isotopomer  the  simulation  was 
based  on  the  same  geometry  and  temperature  parameters 
as  inferred  for  33B3  but  using  a  conventional  Hamiltonian 
of  an  asymmetric  top  molecule.  It  was  also  assumed  that 
only  the  B2-Ax  transition,  which  is  a  type,  occurs.  The 
rotational  constants  used  for  the  32B3  molecule  were 
A:B:C=  1.2696: 1.1906:0. 6144cm- 1  for  the  ground,  and 
1.2457:1.1682:0.6029  cm~'  for  the  excited  state.  The  spec¬ 
trum  of  32B3  was  then  shifted  by  —1.94  cm-1  (from  7  00  of 
33B3)  to  get  the  best  agreement  with  the  observation  (Fig.  2). 
The  isotopic  shift  of  —  1 .94  cm~ 1  is  reasonable  for  an  origin 
band. 

In  all,  34  rotational  transitions  were  assigned;  14  of  them 
are  attributed  to  the  33B3  molecule  and  the  rest  to  32B3  (Table 
III).  The  presented  rotational  assignment  and  molecular 
structure  of  the  B3  molecule  should  be  considered  as  tenta¬ 
tive,  because  the  insufficient  S/N  ratio  of  the  spectrum  de¬ 
creases  the  accuracy  of  the  molecular  structure  determina- 


Cias  et  at. 


TABLE  III.  Frequencies  (in  cm  *)  and  assignments  in  the  2  2E'<—  X2A[ 
absoiption  spectrum  of  32B3  .“ 


NKaKc 

^obs 

^calc 

^obs  ^calc 

2 12" 3 13 

21  847.174 

21  847.197 

-0.023 

2 02" 3 03 

21  847.183 

-0.009 

3  13_4  14 

21  845.859 

21  845.877 

-0.018 

^ 03" ^ 04 

21  845.876 

-0.017 

4[4‘5i5 

4  04"  5  05 

21  844.525 

21  844.544 

-0.019 

5l5"6i6 

5o5"6o6 

21  843.220 

21  843.188 

0.032 

6i6‘7  n 

6  06"  7  07 

21  841.821 

21  841.809 

0.012 

7i7-8is 

7  07" 8  os 

21  840.420 

21  840.407 

0.013 

4 i4‘ 3 13 

21  856.843 

21  856.832 

0.011 

4 04" 3 03 

21  856.833 

0.010 

5l5‘4l4 

5o5‘4o4 

21  857.952 

21  857.934 

0.018 

6l6"5i5 

6 06" 5 05 

21  859.041 

21  859.013 

0.028 

7  17‘6l6 
7o7‘6o6 

21  860.125 

21  860.068 

0.057 

“The  assignment  and  calculated  frequencies  were  derived  on  the  assumption 
that  the  bond  lengths  of  32B3  and  33B3  are  equal,  and  the  isotopic  shift  is 
-1.94  cm-1. 


tion.  Lines  of  3IB3  (=H9%  of  the  33B3  intensity)  disturb  the 
line  profiles  of  32B3  and  33B3  transitions.  There  are  few  un¬ 
assigned  lines,  few  missing  ones,  and  some  are  obscured  by 
intense  neon  atomic  lines  (Fig.  2).  However,  the  consistency 
of  the  whole  analysis  is  sufficient  to  conclude  that  the 
observed  rotationally  resolved  spectrum  is  indeed  the 
2  2 7s  X  2 A  J  transition  of  the  cyclic  B3  molecule,  and  leads 
to  the  first  structural  information  in  the  gas  phase. 
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Abstract 

The  1 2A"  <—  X2A'j  electronic  transition  of  the  B3  molecule  with  origin  band  at  5990  cm-1  has  been  observed  in  a  6  K  neon 
matrix.  A  vibrational  progression  in  the  spectrum  corresponds  to  the  excitation  of  the  v^a',)  vibrational  mode  in  the  1  2A"  state 
with  a  frequency  of  ~1092  cm'* 1.  The  band  system  was  detected  after  laser  ablation  of  a  boron  rod  and  the  origin  band  also  after 
mass-selected  deposition  of  BJ  anions  with  subsequent  neutralization. 

©  2005  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Structure  and  properties  of  boron  clusters  are  inter¬ 
esting  for  fundamental  science.  Boron  is  the  only  ele¬ 
ment  except  carbon  that  can  build  molecules  of  any 
size  by  covalent  bonds,  and  there  are  a  number  of  appli¬ 
cations  [1,2],  Nevertheless,  even  molecules  as  small  as  B3 
are  spectroscopically  not  thoroughly  investigated.  The 
first  study  of  B3  was  by  ESR  [3]  followed  by  IR  spectros¬ 
copy  [4]  both  in  rare  gas  matrices.  A  previous  publica¬ 
tion  reported  the  12E'  <—  X2A|  and  22E'  <—  X2A'j 
electronic  transitions  of  B3  in  a  6  K  neon  matrix  [5]. 
The  rotationally  resolved  origin  band  of  the 
22E'  <—  X 2  A,  transition  has  been  measured  since  then 
in  absorption  in  the  gas  phase  and  analyzed  [6].  Ab  ini¬ 
tio  calculations  on  the  ground  and  excited  electronic 
states  of  B3  are  available  [5,7].  In  the  present  work  the 
1 2A"  4—  X2Aj  electronic  transition  of  B3  has  been  ob¬ 
served  for  the  first  time  in  a  6  K  neon  matrix. 
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2.  Experimental 

The  main  features  of  the  experimental  set  up  have 
been  described  [8].  Two  approaches  were  used  to  pro¬ 
duce  boron  clusters:  one  with  a  cesium  sputter  ion 
source  and  the  other  by  laser  ablation.  Boron  rod 
(99.6%  purity)  with  a  natural  isotopic  composition  was 
used  as  a  target  for  the  laser  vaporization.  B4C,  LaB6 
and  BN  rods  were  used  as  well. 

A  slowly  rotating  sample  rod  was  installed  ~10cm 
from  the  matrix  substrate.  A  532-nm  pulsed  Yag  laser 
beam,  ~20  mJ,  20  EIz,  was  focused  on  the  rod.  A  mix¬ 
ture  of  ionic  and  neutral  ablation  products  was  co¬ 
deposited  with  excess  of  neon  (or  argon)  onto  a 
sapphire  substrate  coated  with  rhodium  and  cooled 
to  a  temperature  of  6  K.  Spectra  in  the  11  GO- 
12  000  cm'1  range  were  recorded  after  sampling  the 
neon  matrix  by  a  Fourier  transform  spectrometer 
applying  a  double  reflection  technique.  The  light  beam 
reflects  from  the  metal  surface  of  the  substrate  on 
which  the  neon  matrix  sits,  then  from  a  mirror  into 
the  matrix  again  and  back  from  the  substrate. 
Absorption  spectra  in  the  220-1100  nm  region  have 
been  measured  by  passing  monochromatic  light 
through  the  matrix  in  a  wave-guide  manner  [9]. 
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In  order  to  assign  the  observed  band  system,  experi¬ 
ments  with  mass-selected  boron  anions,  produced  in  a 
cesium  sputter  source  [10]  have  been  carried  out.  The 
anions  produced  in  the  source  were  guided  by  means 
of  electrostatic  lenses  into  the  quadrupole  mass  filter 
and  were  separated  from  neutral  molecules  by  deflecting 
the  beam  through  90°.  The  Bf  ions  were  mass  selected 
and  co-deposited  with  excess  of  neon  to  form  a  matrix 
at  ~6  K.  A  typical  current  of  Bf  impinging  on  the  ma¬ 
trix  was  ~5  nA.  This  technique  allowed  the  deposition 
of  BJ  exclusively  but  it  was  not  possible  to  obtain  a  high 
concentration  of  neutral  B3  in  the  matrix. 

Neutral  B3  was  produced  by  irradiation  of  the  matrix 
with  a  mercury  medium  pressure  lamp  during  ~30  min. 
An  important  procedure  was  annealing  up  to  8.5  K 
making  the  neon  matrix  softer  and  thus  allowing  relax¬ 
ation  of  the  trapped  molecules  to  their  energetic  minima 
(and  even  some  reactions  to  proceed). 

3.  Observations 

Several  new  bands  were  seen  after  laser  vaporization 
of  the  boron  rod  in  the  infrared  (Fig.  1),  in  addition  to 
the  previously  reported  12E'  +—  X2A]  and  22E' <— 
X2A',  absorption  systems  of  the  B3  molecule  [5,6]  and 
those  of  the  B  and  B2  species  [11-13]  in  the  visible  and 
UV.  The  new  system  comprises  three  peaks,  separated 
by  ~1092  cm-1,  with  origin  band  at  5990  cm-1  (Table 
1).  This  was  seen  after  vaporization  of  B4C  and  LaB6 
rods  as  well.  The  intensity  of  the  bands  increases  after 
irradiation  of  the  matrix  with  UV  photons  indicating  a 
neutralization  process.  Thus  one  can  conclude  that  the 
carrier  is  a  neutral  boron  molecule.  The  band  system 
is  not  likely  due  to  a  boron  and  hydrogen  (common 


Fig.  1.  The  12A('  <—  X2A]  electronic  absorption  spectrum  of  B3  in  a 
6  K  neon  matrix  after  laser  vaporization  of  a  pure  boron  rod.  The 
band  marked  with  *  is  due  to  an  unknown  impurity. 


Table  1 

Positions  of  the  band  maxima  (±0.2  cm-1)  in  the  12A"  <—  X2A] 
electronic  transition  of  B3  in  a  6  K  neon  matrix 


v  (cm  ) 

Av  (cm  *) 

Vibrational  transition 

5990 

0 

og 

7085 

1095 

lj 

8169 

2179 

lo 

Fig.  2.  The  origin  band  of  the  1 2  Aj  <—  X 2  A]  electronic  absorption  of 
B3  in  a  6  K  neon  matrix  after  laser  vaporization  of  a  pure  boron  rod 
(lower  trace),  and  after  mass-selected  deposition  of  B3  with  subsequent 
neutralization  by  UV  irradiation  (upper  trace  -  multiplied  by  a  factor 
of  50  and  smoothed). 


impurity)  containing  species.  When  1%  of  hydrogen 
was  introduced  into  the  matrix  (Ne/H2  ratio  of  100/1) 
the  intensity  of  bands  remained  the  same.  The  absorp¬ 
tion  band  at  6268  cm- 1  (marked  with  *)  is  due  to  an  un¬ 
known  impurity.  Its  intensity  depended  on  the  rod  used. 

Fine  structure  of  the  bands  is  due  to  matrix  site  ef¬ 
fects.  The  trapped  molecules  have  different  orientations 
in  the  neon  lattice.  One  can  see  that  the  three  bands  of 
the  new  system  possess  the  same  structure.  After  anneal¬ 
ing  to  8.5  K  intensity  redistribution  is  observed.  The 
site-pattern  is  different  when  argon  is  used  as  the  matrix. 

Mass  selected  deposition  of  B3  followed  by  neutral¬ 
ization  led  also  to  the  appearance  of  the  5990  cm-1 
absorption  band  (Fig.  2).  The  wavelength  of  the  band 
and  the  site  structure  are  exactly  the  same  as  in  the  laser 
vaporization  spectrum.  However,  the  intensity  is  about 
50  times  lower  for  the  mass  selected  case.  The  bands 
at  7085  and  8169  cm-1  corresponding  to  vibrational 
excitation  in  the  upper  electronic  state  were  not  ob¬ 
served  in  the  mass-selected  experiment  due  to  a  low  sig¬ 
nal  to  noise  ratio.  It  is  concluded  that  the  observed 
system  with  the  origin  band  at  5990  cm-1  is  an  elec¬ 
tronic  transition  of  B3. 
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4.  Discussion 

It  is  theoretically  established  that  B3  has  an  equilat¬ 
eral  triangle  structure  (D3h  group)  with  2 A,  symmetry 
in  the  ground  state  with  electronic  configuration 
(la'1)2(le')4(lai[)2(2a,1)1  [5,7,14],  The  only  excited  states 
accessible  from  the  ground  state  by  electric  dipole  tran¬ 
sitions  are  2A"  and  2E'.  The  electronic  transitions 
12E'  X2Aj  and  22E'  X2Aj  of  B3  in  a  neon  ma¬ 
trix  have  been  studied  earlier  [5],  The  ab  initio  calcula¬ 
tions  predict  that  the  transition  to  the  lowest  excited 
state  with  2 Aj  symmetry  lies  at  ~6300  cm- 1 ;  in  an  ear¬ 
lier  calculation  7000  cm-1  was  obtained  [7],  These  are 
energetically  in  good  agreement  with  the  observed  origin 
band  (Fig.  1).  Consequently  the  new  absorption  system 
is  assigned  to  the  12A"  <—  X2Aj  dipole-allowed  elec¬ 
tronic  transition  of  B3. 

Bf  has  been  studied  by  photoelectron  spectroscopy 
[14]  giving  the  vertical  detachment  energies  of  2.82  and 
3.56  eV  to  the  B3  final  states  X 2 A]  and  1 2  A!,'.  The 
0.74  eV  (~6000cm-)  difference  between  these  is  also 
in  agreement  with  the  present  observations. 

Equilateral  triangular  B3  molecule  has  two  normal 
vibrations:  a]  and  e'.  Only  transitions  to  the  totally  sym¬ 
metric  vibrational  levels  of  the  upper  state  ( 1 2  A")  are  al¬ 
lowed  from  the  X2Aj  ground  state,  the  solely  populated 
level  at  6  K.  Hence  the  bands  at  7085  and  8169  cm-1 
(Table  1)  are  assigned  as  the  first  two  members  of  the 
progression  structure  associated  with  the  a]  normal 
mode  excitation  in  the  upper  PA],'  state.  The  V|  (a', )  fre¬ 
quency  is  thus  inferred  to  be  1092  ±  3  cm-  ,  in  agree¬ 


ment  with  the  analysis  of  the  photoelectron  spectrum 
[14]  of  Bf  which  yielded  the  value  1130  ±  30  cm-1. 

Acknowledgments 

This  study  was  supported  by  the  Office  of  Aerospace 
Research  and  Development  and  the  Swiss  National  Sci¬ 
ence  Foundation  (Project  200020-100019/1). 

References 

[1]  R.N.  Grimes,  J.  Chem.  Ed.  81  (2004)  658. 

[2]  L.H.  Hanley,  J.L.  Whitten,  S.L.  Anderson,  J.  Am.  Chem.  Soc.  92 
(1988)  5803. 

[3]  Y.M.  Hamrick.  R.J.  Van  Zee,  W.J.  Weltner,  J.  Chem.  Phys.  96 
(1992)  1767. 

[4]  S.  Li,  R.J.  Van  Zee,  W.J.  Weltner,  Chem.  Phys.  Lett.  262  (1996)  298. 

[5]  M.  Wyss,  E.  Riaplov,  A.  Batalov,  J.P.  Maier,  T.  Weber,  W. 
Meyer,  P.  Rosmus,  J.  Chem.  Phys.  119  (2003)  9703. 

[6]  P.  Cias,  M.  Araki,  A.  Denisov,  J.P.  Maier,  J.  Chem.  Phys.  121 
(2004)  6776. 

[7]  R.  Hernandez,  J.  Simons,  J.  Chem.  Phys.  94  (1991)  2961. 

[8]  J.P.  Maier,  Chem.  Soc.  Rev.  26  (1997)  21. 

[9]  V.E.  Bondybey,  T.J.  Sears,  J.H.  English,  T.A.  Miller,  J.  Chem. 
Phys.  73  (1980)  2063. 

[10]  D.  Forney,  J.  Fulara,  P.  Freivogel,  M.  Jakobi,  D.  Lessen,  J. 
Chem.  Phys.  103  (1995)  48. 

[11]  W.R.M.  Graham,  W.  Weltner,  J.  Chem.  Phys.  65  (1976)  1516. 

[12]  S.  Tam,  M.  Macler,  M.E.  DeRose,  M.E.  Fajardo,  J.  Chem.  Phys. 
113  (2000)  9067. 

[13]  P.J.  Bruna,  J.S.  Wright,  J.  Phys.  B:  At.  Mol.  Opt.  Phys.  23  (1990) 
2197S. 

[14]  H.-J.  Zhai,  L.-S.  Wang,  A.N.  Alexandrova,  A.I.  Boldyrev,  V.G. 
Zakrzewski,  J.  Phys.  Chem.  A  107  (2003)  9319. 


HTITlL  ABSTRACT  LITlKS 


THE  JOURNAL  OF  CHEMICAL  PHYSICS  122,  034301  (2005) 

Vibrations  in  the  B 4  rhombic  structure 

Roberto  Linguerri,  Isabelle  Navizet,  and  Pavel  Rosmus 

Laboratoire  de  Chimie  Theorique,  Universite  de  Marne  la  Vallee,  F-77454  Champs  sur  Marne,  France 

Stuart  Carter 

Department  of  Chemistry,  University  of  Reading,  Reading  RG6  2AD,  United  Kingdom 

John  P.  Maier 

Department  of  Chemistry,  University  of  Basel,  Klingelbergstrasse  80  CH-4056  Basel,  Switzerland 

(Received  12  July  2004;  accepted  12  October  2004;  published  online  28  December  2004) 

A  double  minimum  six-dimensional  potential  energy  surface  (PES)  is  determined  in  symmetry 
coordinates  for  the  most  stable  rhombic  (D2h)  B4  isomer  in  its  1Ag  electronic  ground  state  by  fitting 
to  energies  calculated  ab  initio.  The  PES  exhibits  a  barrier  to  the  D4h  square  structure  of  255  cm-1. 
The  vibrational  levels  (7  =  0)  are  calculated  variationally  using  an  approach  which  involves  the 
Watson  kinetic  energy  operator  expressed  in  normal  coordinates.  The  pattern  of  about  65  vibrational 
levels  up  to  1600  cm~'  for  all  stable  isotopomers  is  analyzed.  Analogous  to  the  inversion  in 
ammonia-like  molecules,  the  rhombus  rearrangements  lead  to  splittings  of  the  vibrational  levels.  In 
B4  it  is  the  BXg(D4h)  mode  which  distorts  the  square  molecule  to  its  planar  rhombic  form.  The 
anharmonic  fundamental  vibrational  transitions  of  lB4  are  calculated  to  be  (splittings  in 
parentheses):  G(0)  =  2352(22)  cm-1,  v1  (A  lg)  =  1 136(24)  cm~  ',  ^2(^i^)  =  209(144)  cm-1, 

v3(B2g)  =  1198(19)  cm"1,  v4(fl2u)  =  271(24)  cm"1,  and  v5(Eu)  =  1030(  166)  cnT 1  (D4h 
notation).  Their  variations  in  all  stable  isotopomers  were  investigated.  Due  to  the  presence  of  strong 
anharmonic  resonances  between  the  Blg  in-plane  distortion  and  the  B2u  out-of-plane  bending 
modes,  the  higher  overtones  and  combination  levels  are  difficult  to  assign  unequivocally.  ©  2005 
American  Institute  of  Physics.  [DOI:  10.1063/1.1828045] 


I.  INTRODUCTION 

Boron  catenation  leads  to  a  large  variety  of  stable  struc¬ 
tures  which  are  of  interest  in  many  research  areas  and 
applications.1  They  have  been  the  subject  of  numerous 
experimental  and  theoretical  investigations. 

Spectral  data  are  still  limited  for  pure  boron  clusters. 
Two  electronic  transitions  of  the  B2  molecule  have  been  ob¬ 
served  in  the  gas  phase8  and  in  cold  matrices.9  Recently,  the 
electronic  absorption  spectra  of  fi3  and  Bf  in  Ne  matrix,29 
and  the  gas  phase  photoelectron  spectra30  of  Bf  and  Bf 
were  reported. 

The  present  work  deals  with  the  ab  initio  calculations  of 
the  low  lying  vibrational  levels  in  the  most  stable  rhombic 
B4  isomer  in  its  X  1 A ,,  ground  state.  In  previous  theoretical 
studies10-12’22^28’30  the  equilibrium  structures11’12’22-25’28’30 

9<r  in 

and  harmonic  wave  numbers  ’"  ’  were  calculated  for  the 
Dlh  rhombic  minimum  or  for  the  square  D4h  transition  state. 
However,  the  computations  of  a  realistic  pattern  of  vibra¬ 
tional  levels  in  B4  require  methods  going  beyond  the  har¬ 
monic  approximation.  The  existence  of  two  equivalent 
minima  (there  are  six  equivalent  ones  if  permutations  includ¬ 
ing  bond  breaking  are  considered  to  be  feasible  symmetry 
operations)  leads  to  splittings  and  additional  complications 
for  the  vibrational  assignments  are  caused  by  the  presence  of 
strong  low  lying  anharmonic  resonances. 

In  the  tetraatomic  cyclic  clusters,  the  lower  symmetry 
minima  can  correlate  with  electronically  degenerate  states  of 
higher  symmetry  forming  Jahn-Teller  pairs.  In  fact,  the 
theory  for  this  type  of  Born-Oppenheimer  breakdown  has 
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been  originally  formulated  for  the  X4  planar  molecules.31,32 
Based  on  the  model  of  Pearson,33  the  planar  rhombic  equi¬ 
librium  structures  have  been  considered  in  previous  theoret¬ 
ical  studies  (for  instance,  in  Ref.  25)  as  a  result  of  the 
second-order  (pseudo)  Jahn-Teller  effect.  The  coupling  of  at 
least  two  electronic  states  is  supposed  to  lead  to  the  D2h 
rhomboid  minimum.  The  archetypical  example33  of  the 
pseudo-Jahn-Teller  effect  is  the  ammonia-like  D3h  / C3v 
structure,  in  which  the  ground  and  the  first  excited  singlet 
states  are  coupled  by  the  inversion  mode.  Recent  calculations 
clearly  showed  that  the  near-equilibrium  vibrational  problem 
can  be  solved  very  accurately  within  the  Bom-Oppenheimer 
approximation.34  36  In  B4  the  coupling  is  via  the  B  lg  mode, 
which  distorts  the  square  molecule  to  its  planar  rhombic 
X  lA  „  form.  Like  in  ammonia,  the  first  electronically  excited 
singlet  state  in  B4  is  energetically  well  separated  from  the 
ground  state.  We  have  calculated  the  vertical  excitation  en¬ 
ergy  to  the  A  1 Z? !  „  state  at  3.1  eV  above  the  rhombic  equi¬ 
librium  and  3.3  eV  above  the  square  barrier  (cf.  following 
section).  In  the  photodetachment  electron  spectra  of  Bf  (Ref. 
30)  only  excited  triplet  states  were  detected,  lying  between  1 
and  2.4  eV  above  the  X  ]  A  g  ground  state.  The  excited  singlet 
states  were  calculated  at  higher  energies.30  On  the  electronic 
ground  state  potential  energy  surface  (PES)  the  next  higher 
isomer  having  C2v  symmetry  has  been  calculated  to  he  at 
1.77  eV,11  or  2.33  eV12  with  a  barrier  on  the  isomerization 
path  of  0.43  eV12  above  the  rhombic  structure.  Therefore,  for 
the  energy  region  of  the  fundamental  vibrational  transitions  a 
single  adiabatic  PES  with  two  equivalent  minima  is  consid- 
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TABLE  I.  Definition  of  the  symmetry  relations  of  the  planar'  structures  used 
in  the  assignments. 


Mode 

^4  h 

10n  1 1  D 

Ij  4,  D  4 

£2/ 

\0d  1 1  d 

O  2  09 

C2„b 

10bus3,  ub10b3 

c  b 

H)d  11d 

D  2  O  2 

Cs 

1 

Af 

^1 

^1 

A' 

2 

B>s 

^1 

B  2 

A' 

3 

S2s 

B  2 

^  I 

A’ 

4 

Biu 

Bin 

B 1 

a2 

A" 

5 

Eu 

®2u  +  ^3u 

a,+b2 

a1+b7 

2  A' 

az  axis  is 
bC2  lies 

perpendicular  to  the  (r(xy)  molecular  plane, 
in  the  z  axis  and  <r(yz)  is  the  molecular  plane. 

ered  to  be  a  sufficiently  good  approximation. 

Similar  rhombus/square  PES  topology  has  been  pre¬ 
dicted  in  electronic  ground  and  excited  states  of  numerous 
X4  clusters,  for  instance,  in  Li4, 37,38  C4,39  or  Au4. 40,41  Inter¬ 
estingly,  in  Si^C  two  distinct  rhombic  planar  structures  were 
detected  experimentally42  and  one  of  the  isomers  in  interstel¬ 
lar  space.42 

The  nuclear  motion  problem  has  been  solved  by  the 
MULTIMODE  approach.43  It  has  been  demonstrated  that  MUL¬ 
TIMODE  can  deliver  a  near-exact  set  of  vibrational  energies 
for  molecules  such  as  ammonia,44  where  tunneling  through  a 
potential  barrier  is  experienced.  Using  the  same  PES,  identi¬ 
cal  variational  results  were  obtained  with  MULTIMODE  with 
the  full  Hamiltonian  in  internal  coordinates45  and  with  the 
reaction  path  Hamiltonian.46  The  MULTIMODE  method  re¬ 
quires  a  stationary  point  at  which  a  normal  mode  analysis  is 
performed.  The  obvious  choice  for  ammonia,  where  tunnel¬ 
ing  splittings  were  of  interest,  was  the  inversion  saddle  point 
of  Dih  symmetry.  Similarly,  for  B4 ,  the  square  D4h  saddle 
point  is  used  to  obtain  the  normal  modes  required  for  the 
variational  calculations  of  the  anharmonic  wave  functions. 

II.  COMPUTATIONAL  APPROACH 

All  electronic  structure  calculations  were  carried  out  us¬ 
ing  the  MOLPRO  program.47  The  valence  electronic  configu¬ 
ration  of  the  B4  electronic  ground  state  X  *,4  „  ( /J 2 y, )  is 
la2g\b\u\b\u\b\u\b\g2a2g.  Because  it  dominates  the  con¬ 
figuration  expansion  in  the  region  near  the  D4h  barrier  and  at 
the  minimum  Dlh  structure,  we  used  the  coupled  cluster  with 
perturbative  triples  (CCSD(T))  method48  and  the  correlation- 
consistent  polarized  valence  quadruple  zeta  basis  set 
(cc-pVQZ)  of  Dunning49  to  calculate  the  643  energies  used 
in  the  analytical  fit.  The  core  and  core-valence  electron  cor¬ 
relation  was  neglected.  The  vertical  excitation  energy  of  the 
first  electronically  excited  A  1Blg  state  at  the  optimized 
rhombic  equilibrium  structure  (cf.  following  section)  has 
been  calculated  to  be  3.12  eV  equation  of  motion-coupled 
clusters  (EOM-CC);50  cc-pVQZ  basis  set49  and  2.92  eV  full 
valence  complete  active  space  self-consistent  field 
(CASSCF),51  aug-cc-pVTZ  basis  set49,  for  the  square  struc¬ 
ture  3.35  eV  (EOM-CC)  and  3.82  eV  (full  valence 
CASSCF).  As  we  have  performed  calculations  for  all  six 
stable  isotopomers  of  B4  we  define  for  clarity  the  symmetry 
relations  between  different  point  groups  used.  In  Table  I  the 
definition  of  axis,  planes,  and  corresponding  irreducible  rep¬ 
resentations  are  given.  After  having  tested  various  types  of 


FIG.  1.  Nuclear  distances  employed  in  the  definition  of  the  symmetry  coor¬ 
dinates. 


coordinates  the  best  six-dimensional  fit  was  obtained  with 
the  following  D2h  symmetry-adapted  coordinates  (see  Fig.  1 
for  the  definition  of  the  distances): 

5;  =  4(  AR  |  +  AR2  +  AR2  +  AR4), 

Ag 

S2  =  ^(ARl  — AR2  + AR2  — AR4), 

Bis 


S3  =^(AR1-ARt-AR2  +  AR4), 

B2u 


S4  — ^(  AR  [  +  ARo  —  AR2  —  A/?4), 

B3u 


S5  =^(AR5-AR6), 

A*  V2 


=  At, 


t being  the  dihedral  angle  between  the  planes  Bl-B1-B 3  and 
B2-B2-B4.  The  following  polynomial  expansion  is  used  to 
represent  the  six-dimensional  double  minimum  PES: 

4 


V(SUS2,...,S6)=^  2  cIJX,3mt„s\sj2sk3slAs^sn6. 


i,j,k,l,m,n  =  0 


In  the  expansion  only  the  totally  symmetric  products  with 
i  +j  +  k  + 1  +  m  sS4  were  retained.  In  the  planar  square  refer¬ 
ence  geometry  /?°  =  2.8986  bohrs  (i  =  1,4),  /??  =  2.8986 

X  ^2  bohr(y  =  5,6),  and  r°  =  0.  The  following  relation  in¬ 
cluding  permutations  of  the  coordinates  holds: 
V(Sl,S2,S3,S4,S5,S6)  =  V(Sl,S2,-S4,-S2,-S5,-S6). 
The  goal  was  to  obtain  a  PES  valid  close  to  the  minimum 
and  the  barrier.  The  analytic  form  of  the  double  minimum 
PES  is  available  from  the  data  base.52 

In  Fig.  2  contour  plots  are  displayed  for  some  typical 
cuts  of  the  PES.  The  units  on  the  axis  are  bohr  for  the  dis¬ 
tances  and  degree  for  the  S6  coordinate.  For  a  given  pair  of 
coordinates  all  the  others  are  fixed  at  the  optimized  geometry 
of  the  square  transition  state  (cf.  Table  II).  The  two- 
dimensional  cut  of  the  S4  versus  S2  coordinates  shows  an 
isotropic  form.  In  fact,  due  to  the  totally  symmetric  products 
in  the  expansion  and  the  aforementioned  symmetry  relation, 
the  cut  must  be  symmetric  with  respect  to  S2  =  —S4  and  A3 
=  S4  coordinate  permutations.  The  second  cut,  S(<  versus  S , , 
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FIG.  2.  Four  contour  plots  of  the  double  minimum  six-dimensional  PES  of 
the  S4  rhombic  cluster  in  symmetry  coordinates.  The  distance  is  in  bohr  and 
angle  in  degree.  Contour  a:  energy  range  6000  cm~\  contour  step  250  cm-1. 
Contour  b:  energy  range  4000  cm-1,  contour  step  250  cm-1.  Contour  c: 
energy  range  2000  cm-1,  contour  step  50  cm-1.  Contour  d:  energy  range 
1750  cm-1,  contour  step  250  cm-1. 


displays  the  potential  for  all  in-plane  (S6  =  0)  and  out-of- 
plane  (S6¥=0)  geometries,  in  which  Rl=R2  =  Ri  =  R4  and 
R5  =  R6.  The  S 5  versus  Sl  cut,  including  the  rhombic 
minima,  gives  information  on  how  the  equal  sides  of  the 
tetragon  (S2  =  S2  =  S4  =  0)  become  larger  when  passing  from 
the  square  to  the  rhombic  structure  (the  S  j  value  is  positive 
for  the  rhombic  minimum  and  zero  for  the  square  structure). 
The  S 6  versus  S 5  cut  in  the  last  diagram  shows  the  change  of 
the  energy  barrier  with  increasing  out-of-plane  angle,  for 
squeezing  a  tetragon  with  Rl=R2  =  R3  =  R4  =  R®. 

The  barrier  height  (on  the  adiabatic  PES)  relative  to  the 
rhombic  minimum  has  been  calculated  to  be  255  cm-1.  Pre¬ 
viously  Martin  et  al.25  reported  178  cm~'  (G1  corrected 
value),  Zhai  et  al.30  2  8  3  cm~'  (B3LYP),  and  Kato  and 
Tanaka26  1014  cnT1  (MP4). 


The  analysis  of  the  PES  indicates  that  the  path  corre¬ 
sponding  to  two  clockwise  and  two  anticlockwise  in-phase 
nuclear  in-plane  rotations  for  the  r  cos(  6)  coordinate  with  r 
=  (dl—d2)/4,  and  d\  being  the  long  and  d2  the  short  diago¬ 
nal  of  the  rhombus,  represents  a  good  approximation  for  the 
periodic  Blg  ring  distortion.  However,  there  is  no  path  con¬ 
necting  the  two  rhombic  minima  below  the  square  transition 
state. 

The  PES  described  above  was  used  in  vibrational  varia¬ 
tional  calculations  with  the  code  MULTIMODE.53  This  proce¬ 
dure  utilizes  the  exact  Watson  kinetic  energy  operator  in  nor¬ 
mal  coordinates,54  and  therefore  some  suitable  reference 
structure  at  which  the  3N—6  normal  vibrations  are  defined 
must  be  first  located.  In  order  to  be  even  handed  about  the 
two  equivalent  minima  on  the  PES,  the  choice  for  such  a 
reference  structure  is  that  of  the  square-planar  transition 
state.  The  normal  mode  vectors  at  the  transition  state  can 
then  be  readily  specified  in  terms  of  the  irreducible  represen¬ 
tations  of  the  symmetry  point  group  for  any  chosen  isoto- 
pomer.  However,  due  to  the  specific  difficulties  in  employing 
C2v  symmetry  (the  highest  currently  available  in  MULTI- 
MODE)  for  molecules  with  E  modes  in  the  D4h  point  group, 
all  we  can  do  is  to  use  the  reflection  operation  in  order  to 
symmetry  block  the  configuration  interaction  (Cl)  matrix;  we 
therefore  work  in  Cs  symmetry  throughout  (two  blocks  A' 
and  A"),  and  use  the  vector  symmetries  in  order  to  reclassify 
the  energy  levels  into  irreducible  representations  of  the  point 
group  of  the  transition  state.  The  great  advantage  of  the  MUL¬ 
TIMODE  approach  in  the  application  to  the  B4  problem  is  the 
possibility  to  designate  the  anharmonic  vibrational  states  in 
terms  of  linear  combinations  of  vibrational  configurations 
with  harmonic  quantum  numbers.  This  is  a  universal  bonus 
gained  from  using  normal  coordinates  and  has  been  the 
method  of  choice  for  spectroscopists  for  many  years. 

The  methodology  adopted  in  the  MULTIMODE  approach 
has  been  fully  described  on  a  number  of  occasions  (see  Ref. 
53  and  references  therein).  We  therefore  only  outline  the  ba¬ 
sic  features  relevant  to  this  work.  The  initial  step  in  the  pro¬ 
cedure  is  to  optimize  a  set  of  36  harmonic-oscillator  primi- 


TABLE  II.  Optimized  geometries  (diagonal  distances)  and  the  harmonic  wave  numbers  of  the  fundamental 
vibrational  transition  in  the  electronic  ground  states  of  B4 . 


B 4  rhombus  (Z)2/,) 

B4  square  (D4h) 

This  work3 

b 

c 

This  work3 
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Diagonals 
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1112 
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1136 

(cm-1) 

co2 
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1240 
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Bu 
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w5 

B4„ 

1013 

970 

966 

Eu 

1148 

1138 

"5 

Bt,u 

1235 

1251 

1201 

E„ 

1148 

1138 

aThe  harmonic  wave  numbers  were  calculated  with  the  relative  atomic  mass  of  10.810  for  boron;  cc-pVQZ 
basis  was  used  (Ref.  49),  the  CCSD(T)  total  energies  at  the  rhombic  and  square  structures  were  calculated  to 
be  —98.908  391  a.u.  and  —98.907  260  a.u.  The  standard  molpro  convergency  criteria  for  the  gradients  and 
energies  were  used. 

bMP2  approach/6-3 11G*;  no  relative  mass  of  B  was  given  (Ref.  25). 
cRCCSD(T)/6-311+G*,  Emin=  — 98.839  444  a.u.;  no  relative  mass  of  B  was  given  (Ref.  30). 
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tive  basis  functions  (NBF)  for  3N  —  1  normal  modes  by 
using  effective  one-dimensional  cuts  in  the  PES  relative  to 
the  “squeezing”  mode.  These  one-dimensional  potentials  are 
then  used  in  variational  calculations  with  the  equivalent  one¬ 
dimensional  Hamiltonians,  respectively,  for  which  ten  con¬ 
verged  subsets  (NVF)  are  retained  for  subsequent  use.  The 
raw  one-dimensional,  double -minimum  cut  is  used  for  the 
squeezing  mode  itself  (NBF=46,  NVF=15),  in  order  to  op¬ 
timize  the  corresponding  harmonic-oscillator  basis.  The 
primitive  harmonic-oscillator  basis  is  defined  by  the  values 
of  oj  arising  from  the  Wilson  FG  normal  mode  analysis,  with 
the  exception  of  that  for  the  squeezing  mode,  whose  eigen¬ 
value  is  imaginary  at  the  planar  transition  state.  In  this  case 
pilot  calculations  are  carried  out  with  various  values  of  u> 
assigned  to  the  harmonic -oscillator  basis,  until  an  optimum 
value  is  obtained  which  spans  sufficiently  into  the  high- 
energy  regions  of  the  one-dimensional  cut.  Providing  that  a 
sufficient  number  of  primitive  harmonic-oscillator  functions 
are  used,  fully  converged  one-dimensional  expansion  func¬ 
tions  can  be  obtained  for  modes  with  such  gross  anharmo- 
nicity.  This  optimization  procedure  is  carried  out  using 
Gauss-hermite  numerical  quadrature  involving  quadrature 
points  to  integrate  all  matrix  elements  (50  for  the  3N—1 
normal  modes  and  60  for  the  squeezing  mode),  and  the  re¬ 
sulting  contracted  expansion  functions  are  then  used  to  re¬ 
duce  the  integration  points  via  the  HEG  algorithm.55 

The  optimized  one-dimensional  expansion  set  is  now 
used  in  a  vibrational  self-consistent-field  procedure  for  the 
zero-point  level.  Integration  of  the  matrix  elements  takes 
place  in  a  series  of  steps  (^5),  during  which  the  overcount¬ 
ing  of  each  M-mode  potential  resulting  from  M  + 1  -mode 
cuts  (and  so  on)  is  taken  into  account.43  Fikewise,  integration 
of  the  Coriolis  terms  which  involve  the  inverse  moment-of- 
inertia  tensor  jj,  is  performed  in  an  analogous  way.  For  each 
M-mode  cut  in  the  PES  and  /jl ,  the  equivalent  M-mode 
Hamiltonian  is  adopted.  The  result  of  this  analysis  is  an  im¬ 
proved  one -dimensional  expansion  set  in  the  NVF  basis  for 
each  of  the  3N  —  6  modes,  and  these  are  now  used  in  a  final 
Cl  calculation  to  determine  the  vibrational  energies.  The  Cl 
matrix  itself  is  constructed  by  imposing  a  maximum  quan¬ 
tum  for  each  mode,  together  with  a  maximum  sum-over- 
quanta  for  all  modes,  in  a  series  of  M-mode  blocks  (M 
ss5).  The  complete  Cl  matrix  is  the  total  of  all  such 
blocks.53  Finally,  the  resulting  Cl  matrix  is  split  into  A '  and 
A"  symmetry  blocks.  As  with  all  variational  calculations,  in¬ 
creased  Cl  matrices  are  diagonalized  until  the  required  con¬ 
vergence  is  reached.  In  MULTIMODE  calculations,  a  further 
convergence  criterion  is  required,  and  that  is  for  the  M-mode 
integration  procedure.  In  this  study,  virtually  total  conver¬ 
gence  in  all  levels  was  achieved  for  M  =  5,  and  this  therefore 
removes  the  inherent  approximation  from  the  calculations. 


III.  RESULTS  AND  DISCUSSION 

The  computational  approach  used  for  the  generation  of 
the  six-dimensional  PES  yielded  the  following  spectroscopic 
constants  for  the  nB2X3’Z  state:  Rf=1.593A,  u>e 

_  1  '  ’  _  i 

=  1048 cm  ,  and  coexe=  8.92 cm  .  The  corresponding  ex¬ 


perimental  values8  are  Rf=1.590A,  u>e=  1051.3 cm~ *,  and 
a) cxe  =  9. 35  cm-1.  It  is  difficult  to  estimate  the  accuracy  of 
the  fundamental  transitions  evaluated  from  the  six¬ 
dimensional  PES  of  B4,  but  we  expect  to  have  obtained 
similarly  reliable  results.  The  CCSD(T)  geometries  and  har¬ 
monic  wave  numbers  for  the  minimum  and  the  saddle  point 
of  B4  are  compared  with  previous  theoretical  data  in  Table  II. 
The  similarity  of  the  B-B  equilibrium  distances  in  B2  and 
B4  is  related  to  the  covalent  character  of  the  bonds  in  the 
cluster.  For  D4h  structures  only  one  imaginary  wave  number 
is  calculated  for  the  B  \ ,,  mode,  which  transforms  the  square 
transition  state  to  the  rhombic  minimum.  Our  optimized  di¬ 
agonal  distances  and  harmonic  wave  numbers  for  B4  differ 
from  those  calculated  previously.  The  present  results  are  seen 
to  be  more  accurate  due  to  the  better  quality  of  the  electronic 
structure  calculations. 

The  CCSD(T)  values  for  the  cd2  ring  in-plane  distortion 
and  the  u>4  ring  out-of-plane  bending  are  calculated  to  be 
almost  identical.  This  implies  that  overtones  or  combination 
levels  with  even  vibrational  quanta  will  form  anharmonic 
resonances  considerably  complicating  the  splitting  patterns 
and  assignments.  Within  a  given  symmetry  block  other 
nearby  states  interact  as  well. 

In  nature  there  are  two  stable  boron  isotopes,  1  °B  (abun¬ 
dance  19.9%)  and  nB  (80.1%),  which  give  rise  to  six  differ¬ 
ent  isotopomers  of  B4  belonging  to  different  point  groups. 
The  MULTIMODE  calculations  were  all  done  using  Cs  symme¬ 
try  with  the  plane  in  which  the  four  boron  atoms  lie  at  the 
transition  state,  defining  the  symmetry  plane.  In  all  six  iso¬ 
topomers  the  vibrational  states  in  Table  III  are  assigned  with 
the  highest  possible  symmetries  of  the  transition  states.  In 
the  first  column  the  consecutive  numbering  of  all  levels  is 
used  and  the  labels  l  (lower)  and  u  (upper)  indicate  the  split 
pairs.  In  the  second  column  the  harmonic  vibrational  con¬ 
figurations  (with  coefficients  larger  than  0.3  in  the  configu¬ 
ration  expansion)  are  given  for  1  lB4 .  The  mode  numbering  is 
defined  in  Table  I.  The  results  are  presented  in  two  blocks — 
Table  III  for  A'  symmetry  and  Table  IV  for  A"  symmetry.  All 
calculated  levels  for  UB4  up  to  about  1500  cm-1  are  dis¬ 
played  in  Fig.  3,  with  assignments  using  the  irreducible  rep¬ 
resentations  of  the  D4h  point  group.  Each  state  is  designated 
by  specifying  the  dominating  harmonic  configurations  in  the 
total  wave  function  expansion.  The  zero-point  energies  are 
given  relative  to  the  electronic  energy  of  the  transition  state 
and  vary  up  to  122  cm-1  for  different  isotopomers.  The 
vibrational  mode  splittings  were  deduced  from  the 
assignments  in  Tables  III  and  IV.  The  anharmonic  funda¬ 
mental  vibrational  levels  of  11 B4  are  calculated  to  be  (split¬ 
tings  in  parentheses):  G(0)  =  2352(22)  cm-1,  vi(Alg) 
=  1136(24)  cm-1,  v2(Blg)  =209(  144)  cm'1,  v3(B2g) 
=  1 198 (  19)  cm-1,  v4(B2u)  =271(24)  cm-1,  and  v5(Eu) 
=  1030(  166)  cm~ 1 .  The  corresponding  values  for  the  second 
abundant  l0BnB3  isotopomer  are  2364  (23),  1148  (24),  212 
(148),  1218  (21),  274  (25),  and  1041,  1042  (158,160)  (all 
values  in  cm-1).  All  IR  or  Raman  allowed  fundamentals  of 
these  two  isotopomers  will  be  clearly  discernible  in  their 
spectra.  According  to  the  splittings  in  the  ground  state  of 
about  22-26  cm"  the  rhombus  will  rearrange  within  about 
0.75  ps,  and  at  room  temperature,  because  the  Blg  and 
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TABLE  III.  Vibrational  pattern  in  the  B4  isotopomers  (A'  symmetry). 


Abundance  (%) 

llB 4 

41.2 

wBnB} 

40.9 

wB2nB2 

10.2 

lOp  Up 

D  2  D  2 

5.1 

l0B3nB 

2.5 

wB4 

0.2 

Assignment3 

T>4  * 

c2„ 

C2v 

E>ih 

c2v 

E>4h 

o' 

G(0) 

4  ig 

2332.76 

A, 

2363.60 

A  i 

2394.13 

Ag 

2394.36 

^1 

2424.58 

2454.73 

0" 

21 

Big 

21.89 

^1 

22.97 

Bi 

24.03 

Ag 

24.12 

^1 

25.18 

Big 

26.32 

1' 

22 

209.15 

^1 

212.07 

A  i 

215.12 

Ag 

215.02 

^1 

218.21 

Ais 

221.46 

1" 

23 

Blg 

353.36 

A, 

359.70 

Bi 

366.30 

Ag 

366.04 

A, 

372.91 

Big 

379.81 

21 

42  +  24 

Alg 

534.62 

A, 

541.69 

^1 

548.57 

Ag 

548.48 

^1 

555.21 

A  ig 

561.76 

21 

24  +  42 

A  i, 

556.36 

A, 

565.28 

^1 

574.68 

Ag 

574.33 

^1 

584.23 

Ai8 

594.27 

2" 

2 1 42 

Big 

566.79 

^1 

574.14 

Bi 

581.37 

Ag 

581.52 

^1 

588.60 

Big 

595.82 

3' 

2242 

724.64 

^1 

733.43 

^1 

742.22 

Ag 

742.18 

^1 

750.99 

759.88 

2“ 

25 

Blg 

765.29 

A, 

778.55 

Bi 

792.32 

Ag 

791.72 

^1 

806.02 

Big 

820.25 

3" 

2342 

Blg 

860.94 

A, 

872.57 

Bi 

884.35 

Ag 

884.17 

^1 

896.09 

Big 

908.12 

4' 

26  +  2442 

988.01 

^1 

1002.77 

^1 

1017.65 

Ag 

1017.13 

^1 

1032.12 

A  ig 

1046.72 

5' 

51  +  2I51 

F-u 

1030.05 

^1 

1041.37 

^1 

1053.05 

b2u 

1054.31 

^1 

1066.36 

Eu 

1080.72 

5' 

51  +  2151 

Eu 

1030.11 

Bi 

1042.35 

Bi 

1054.11 

Blu 

1057.15 

Bi 

1068.28 

Eu 

1080.81 

4' 

2442  +  26  +  44 

Aig 

1036.31 

Ai 

1052.02 

A  i 

1068.21 

AS 

1067.52 

^1 

1084.36 

A  ig 

1101.17 

4' 

44  +  2442 

A  i, 

1094.36 

Ai 

1107.79 

1121.15 

Ag 

1120.98 

^1 

1134.45 

Aig 

1147.90 

4" 

2 1 44 

Big 

1102.23 

•A  I 

1115.63 

Bi 

1128.67 

Ag 

1128.85 

^1 

1141.72 

Big 

1154.72 

6' 

l1 

1135.76 

^1 

1148.43 

^1 

1161.16 

Ag 

1163.19 

^1 

1175.44 

1190.15 

6" 

21 1 1 

Big 

1159.79 

A, 

1172.51 

Bi 

1184.79 

Ag 

1189.93 

A, 

1202.63 

Big 

1218.94 

5" 

51  +  2151 

E„ 

1196.03 

^1 

1199.80 

A  i 

1213.39 

Biu 

1223.80 

Bi 

1227.26 

Eu 

1252.36 

5" 

2'5 1  +  5 1 

Eu 

1196.08 

Bi 

1214.26 

b2 

1222.48 

B2u 

1225.18 

^1 

1242.74 

Eu 

1252.41 

7' 

31 

Blg 

1198.45 

Bi 

1218.12 

^1 

1242.11 

Big 

1225.75 

Bi 

1249.36 

Big 

1255.58 

4" 

2542  +  27  +  2‘44 

Blg 

1206.48 

^1 

1223.05 

b2 

1239.72 

Ag 

1239.24 

^1 

1256.19 

Big 

1272.89 

7" 

2*3 1 

Alg 

1217.94 

Bi 

1239.47 

b2 

1255.99 

Big 

1250.60 

Bi 

1267.22 

Aig 

1279.17 

8' 

2244 

A  i. 

1239.17 

^1 

1253.93 

^1 

1268.42 

Ag 

1268.36 

^1 

1282.89 

1297.43 

4" 

27  +  2542 

BU 

1282.67 

A, 

1303.99 

b2 

1326.08 

Ag 

1325.64 

^1 

1347.35 

Big 

1372.28 

9‘ 

225 1 

Eu 

1302.35 

Bi 

1316.58 

^1 

1327.45 

B2u 

1333.96 

^1 

1349.56 

Eu 

1367.28 

9' 

225 1 

Eu 

1302.60 

^1 

1317.14 

b2 

1336.11 

Biu 

1336.05 

Bi 

1350.94 

Eu 

1367.62 

10' 

22 1 1 

Aig 

1344.13 

^1 

1363.58 

^1 

1382.20 

Ag 

1377.86 

^1 

1397.03 

A  ig 

1411.39 

8“ 

2344  +  2542 

Big 

1378.39 

^1 

1396.31 

b2 

1414.12 

Ag 

1413.88 

^1 

1431.98 

Big 

1450.08 

11' 

2542 

1405.12 

^1 

1423.42 

^1 

1441.83 

Ag 

1441.19 

^1 

1459.89 

A  ig 

1478.86 

12' 

223 1 

B2g 

1408.77 

Bi 

1426.58 

^1 

1444.26 

Big 

1443.22 

Bi 

1460.98 

Big 

1477.24 

9  « 

235 1 

E„ 

1464.38 

b2 

1480.50 

B2 

1493.83 

Biu 

1502.51 

^1 

1518.89 

Eu 

1541.77 

o« 

235 1 

E„ 

1464.59 

Ai 

1482.65 

^1 

1505.27 

B}u 

1503.46 

Bi 

1521.99 

Eu 

1542.01 

10" 

23 1 1 

Big 

1492.89 

^1 

1516.85 

Bi 

1539.38 

Ag 

1533.56 

A, 

1557.49 

Big 

1574.50 

13' 

23  +  2444 

1529.39 

^1 

1551.44 

A  i 

1573.32 

Ag 

1572.68 

^1 

1594.87 

A  ig 

1616.48 

12" 

233 1 

Alg 

1548.39 

Bi 

1570.28 

b2 

1591.98 

Big 

1589.64 

Bi 

1611.41 

Aig 

1631.25 

13' 

23  +  2642  +  2444 

A  ig 

1567.62 

A  1 

1592.47 

A  i 

1619.08 

Ag 

1618.13 

^1 

1645.58 

A  ig 

1673.27 

14' 

425 1  +  2'425 1 

E„ 

1574.33 

^1 

1591.47 

1609.47 

B2u 

1610.52 

^1 

1629.39 

Eu 

1649.15 

14' 

425 1  +  2*425 1 

Eu 

1574.85 

Bi 

1594.93 

Bi 

1611.28 

Biu 

1614.34 

Bi 

1631.45 

Eu 

1649.69 

“The  first  column  lists  pairs  of  levels  split  by  the  B  lg  (D4h)  motion.  The  /  and  u  labels  designate  lower  and  upper  levels,  respectively.  The  vibrational  quantum 
numbers  of  the  big  distortion  mode  in  the  u  and  /  vibrational  configurations  differ  by  one.  The  pairs  are  numbered  according  to  increasing  wave  numbers. 
When  two  levels  correspond  to  the  same  combination  of  vibrational  configurations,  they  are  attributed  to  the  same  pair.  Second  column  gives  the  leading 
harmonic  vibrational  configurations  of  the  wave  function. 


B2u  fundamentals  will  also  be  populated,  the  cycle  will  un¬ 
dergo  a  complex  motion  including  the  out-of-plane  bending. 

The  term  pattern  is  complicated  by  the  presence  of 
strong  anharmonic  resonances.  In  the  A '  block,  for  instance, 
the  harmonic  states  with  two  quanta  in  the  Big  distortion  and 
two  quanta  in  the  out-of-plane  Blu  form  a  lower  (21  in  Table 
III)  anharmonic  pair  at  535  and  556  cm~'  in  M/i4 ,  their 
upper  counterparts  (2")  are  calculated  to  lie  at  567  and  765 
cm-1.  In  the  A"  block  the  31  pair  lies  at  780  and  823  cm-1, 
the  3"  pair  at  836  and  992  cm-  .  The  assignments  of  the  split 
states  becomes  difficult  going  up  the  ladder  and  the  attribu¬ 
tion  to  the  upper  and  lower  levels  in  Tables  III  and  IV  are 
only  tentative.  Similar  resonances  exist  also  in  other  isoto¬ 
pomers.  The  calculations  yield  essentially  identical  vibra¬ 
tional  energies  for  the  two  possible  C2v  rhombic  isotopomers 


wB2uB2,  but  the  anharmonic  couplings  lead  to  shifts  par¬ 
ticularly  in  the  region  of  the  E  modes. 

IV.  SUMMARY  AND  CONCLUSION 

To  date,  only  the  photodetachment  spectrum  of  rhombic 
B4  has  yielded  experimental  information  on  the  ground  state 
of  B4.  The  detected  band  is  broad  with  no  distinguishable 
vibrational  progression/'0 

In  the  present  work  we  used  a  double  minimum  six¬ 
dimensional  CCSD(T)  PES  for  the  lowest  isomer  of  B4  to 
solve  the  anharmonic  nuclear  motion  problem  by  the  MULTI- 
MODE  approach.  The  calculated  splittings  of  the  fundamental 
levels  (around  20  cm~  )  indicate  that  the  rhombic  structure 
rearranges  on  the  time  scale  of  about  0.75  ps.  The  presence 
of  the  splittings,  anharmonic  resonances,  hot  bands,  and 
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TABLE  IV.  Vibrational  pattern  in  the  B4  isotopomers  (A"  symmetry). 


Abundance  (%) 

"Da 

41.2 

10bub3 

40.9 

10n  Un 
d  2  O  2 

10.2 

1 0  D  1 1  D 

D  2  D  2 

5.1 

2.5 

10b4 

0.2 

Assignment3 

Dai, 

Cm 

Cm 

Dm 

Civ 

Dai, 

1' 

4‘ 

Biu 

271.28 

B, 

274.43 

^2 

277.53 

Du, 

277.53 

S, 

280.58 

Biu 

283.62 

1" 

2 1 4 1 

Ah, 

295.18 

B, 

299.42 

B  i 

303.61 

B  In 

303.72 

S, 

307.84 

Ah, 

312.05 

2l 

224 1 

Dm 

467.28 

B, 

473.17 

A  2 

479.13 

Du, 

479.06 

S, 

485.11 

Dm 

491.22 

2" 

234 1 

Ah, 

607.22 

B  i 

616.13 

B  i 

625.21 

Si. 

625.02 

S, 

634.27 

Ah, 

643.60 

3' 

244*  +  43 

^2  u 

780.37 

B, 

791.07 

A  2 

801.78 

Si. 

801.55 

S, 

812.29 

Dm 

822.90 

3' 

43  +  244‘ 

Blu 

822.86 

B  i 

833.33 

A  2 

843.89 

Si. 

843.76 

S, 

854.47 

Dm 

865.25 

3" 

2 1 43 

Ah, 

835.88 

B  i 

846.24 

B  i 

856.40 

fi,„ 

856.56 

S, 

866.54 

Ah, 

876.64 

4' 

2243 

Dm 

981.52 

B, 

993.24 

^2 

1004.87 

Du, 

1004.77 

S, 

1016.40 

Dm 

1028.04 

3" 

2S4‘ 

Ah, 

992.46 

B  i 

1006.94 

Si 

1021.72 

Si. 

1021.29 

S, 

1036.35 

Am 

1051.39 

4  « 

2343  +  254‘ 

Am 

1115.54 

B, 

1130.13 

Si 

1144.74 

Du, 

1144.47 

S, 

1159.28 

Am 

1174.07 

5' 

2641  +  2443 

D  2u 

1197.32 

B, 

1213.67 

^2 

1230.10 

Bu, 

1229.56 

S, 

1246.27 

Dm 

1262.71 

5' 

2443  +  2641  +  4s 

Din 

1289.36 

B  i 

1307.33 

^2 

1325.39 

B,„ 

1325.11 

S, 

1343.75 

Dm 

1362.26 

6' 

4151  +  214151 

Eg 

1304.69 

B, 

1318.96 

Si 

1334.22 

B2g 

1335.00 

Si 

1351.02 

ES 

1367.67 

6' 

4151  +  2‘4151 

Es 

1304.92 

^2 

1321.52 

^2 

1335.22 

1338.43 

^2 

1352.50 

ES 

1368.00 

5" 

2*45  +  2543 

Ah, 

1361.90 

B, 

1379.49 

s, 

1395.76 

B,„ 

1394.47 

S, 

1411.52 

Ah, 

1426.89 

7' 

45  +  2443 

Di„ 

1367.46 

B, 

1384.22 

^2 

1400.32 

Bl„ 

1400.21 

S, 

1416.46 

Dm 

1432.42 

8' 

4 1 1 1 

Dm 

1404.14 

B, 

1419.99 

^2 

1436.04 

Bu, 

1437.81 

S, 

1453.30 

Dm 

1470.65 

7" 

2741  +  2’45  +  2543 

Ai„ 

1422.01 

Bi 

1441.26 

s, 

1459.93 

B,„ 

1459.21 

S, 

1478.49 

Am 

1496.97 

8" 

2 1 4 1 1 1 

Ai„ 

1430.78 

B, 

1445.73 

s, 

1461.06 

S  i„ 

1467.04 

S, 

1481.89 

Am 

1502.25 

9' 

4 1 3 1 

A  i„ 

1457.21 

A  2 

1467.30 

^2 

1483.37 

A„ 

1490.28 

^2 

1502.05 

A  i„ 

1525.70 

6" 

4151  +  2‘4151 

Eg 

1467.38 

B i 

1484.31 

s, 

1497.54 

B2S 

1501.83 

S, 

1521.77 

ES 

1535.96 

6" 

4151  +  2‘4151 

Eg 

1467.87 

A  2 

1489.35 

^2 

1514.14 

B,s 

1502.88 

^2 

1524.04 

Eg 

1536.47 

gu 

2 1 4 1 3 1 

Du, 

1479.48 

^2 

1506.38 

s, 

1525.63 

A„ 

1517.84 

^2 

1538.05 

Bu 

1552.26 

10' 

2245 

^2  u 

1506.03 

B  i 

1525.82 

^2 

1543.72 

Bu, 

1542.09 

s, 

1561.12 

Dm 

1577.91 

5" 

274*  +  2543 

Ah, 

1521.18 

B, 

1543.74 

s, 

1566.79 

Bu, 

1565.99 

S, 

1589.85 

Am 

1613.24 

11' 

224 1 5 1 

E s 

1562.59 

B, 

1580.29 

s, 

1593.11 

B2g 

1600.07 

S, 

1618.21 

ES 

1639.02 

11' 

224151 

Eg 

1563.18 

^2 

1581.59 

^2 

1602.56 

B2S 

1602.41 

A  2 

1620.06 

Eg 

1639.61 

aSee  footnote  in  Table  III. 
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FIG.  3.  Vibrational  energy  pattern  in  UB4  up  to  1500  cm  .  The  vibrational  states  are  assigned  using  irreducible  representations  of  the  D4h  point  group. 
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complex  rotational  structures  with  different  nuclear  spin  sta¬ 
tistics  will  considerably  complicate  the  assignments  of  the 
gas  phase  IR  and  Raman  spectra  at  298  K.  Some  of  the 
dipole  forbidden  transitions  are  allowed  in  the  less  abundant 
isotopomers  and  will  have  their  own  weak  rotational  transi¬ 
tions.  The  theoretical  approach  developed  in  this  work  is 
generally  applicable  for  tetraatomic  rhombic  clusters  not  ex¬ 
hibiting  the  Jahn-Teller  effect. 
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First-principle  full  valence  complete  active  space  calculations  (CASSCF)  are  reported 
for  the  electronic  states  of  the  D2/,  and  C2„  planar  isomers  of  B4,  B4  and  B4.  This 
approach  is  found  to  reproduce  well  the  experimentally  known  transition  energies  for  the 
two  lowest  excited  triplet  states  of  the  most  stable  rhombic  B4  isomer.  The  pattern  of 
electronic  transitions  up  to  about  4.5  eV  for  B4  singlet  and  B4  doublet  states  calculated 
in  similar  way  should  guide  an  UV  spectroscopic  search. 


1.  Introduction 

Studies  of  small  boron  clusters  are  motivated  by  a  desire  to  understand  the 
manifold  of  chemical  structures,  intramolecular  dynamics  and  reactivity  of 
these  covalently  bound  homonuclear  molecules.  There  are  surprisingly  few 
reports  on  the  spectroscopy  of  such  neutral  or  charged  clusters.  The  boron 
trimer  was  investigated  in  low  temperature  matrices  by  ESR  [1],  IR  [2]  and 
UV  [3]  techniques.  So  far  only  its  2  2E'  «—  X  2 A'  electronic  transition  was 
observed  in  the  gas  phase  [4] .  The  B3  and  B4  clusters  in  their  neutral  and 
anionic  forms  were  studied  by  photoelectron  spectroscopy  [5],  The  first  ion¬ 
ization  energy  [6-8]  and  the  electron  affinity  [5]  of  rhombic  B4  are  known. 
The  mass-spectrometric  studies  [9-1 1]  on  ID  provided  information  about  their 
thermochemical  stabilities  and  fragmentation  paths. 

The  present  contribution  deals  mainly  with  electronically  excited  states 
of  B4  and  B4 .  The  aim  is  to  calculate  the  pattern  of  electronic  transitions  in 
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rhombic  and  planar  C2„  isomers  to  guide  the  spectroscopic  search  for  these 
molecules.  Previous  theoretical  works  have  reported  the  ground  state  equilib¬ 
rium  structures  of  various  isomers,  their  relative  stabilities,  harmonic  frequen¬ 
cies,  nuclear-independent  chemical  shifts,  dipole  polarizability,  hyperpolariz¬ 
ability,  and  thermochemical  data  [12-31].  Recently  we  have  investigated  the 
splittings  of  vibrational  states  in  rhombic  B4  [30].  So  far  the  only  theoretical 
study  concerning  the  electronically  excited  triplet  and  singlet  states  of  B4  was 
performed  in  conjunction  with  the  assignment  of  the  B4  photoelectron  spec¬ 
trum  [5]. 

The  boron  tetramer  represents  a  challenging  problem  from  the  point  of 
view  of  electronic  spectroscopy.  Depending  on  its  formation  different  isomers 
could  coexist  and  their  electronic  states  could  lie  in  the  same  UV  region.  The 
states  in  rhombic  isomers  can  correlate  with  electronically  degenerate  states 
of  D4/l  symmetry  forming  Jahn-Teller  pairs.  Perturbations  can  arise  from  vi- 
bronic  couplings  and  spin  multiplet  crossings.  As  in  the  B4  rhombic  ground 
state  [30],  the  minimum/square  transition  state  topology  of  the  potential  en¬ 
ergy  surface  (PES)  can  lead  to  splittings  of  vibrational  modes  complicating  the 
vibronic  coupling  effects. 

The  present  study  shows  that  the  ground  states  of  B4  and  B4  consist  of  two 
vibronically  coupled  electronic  states.  The  most  intense  electronic  transitions 
in  B4  and  B4  will  be  those  to  Jahn-Teller  states.  It  is  found  that  the  pattern  of 
dipole  allowed  electronic  transitions  in  rhombic  and  C2„  isomers  should  allow 
them  to  be  recognized  by  their  UV  spectra. 


2.  Computational  approach 

All  electronic  structure  calculations  were  carried  out  using  the  MOLPRO  pro¬ 
gram  [32].  The  valence  electronic  configurations  for  the  lowest  isomers  of  the 
B4  electronic  ground  state  are: 

X  'A,  (D2,  rhombus)  la2  lb:;,,  lb23u  lb2lu  lb2,,  2a2 
X  ‘A!  (C2„  planar)  la2  2a2  lb2  3a2  lb,  4a2 
X  'E+  (D^u)  lcr2  lcr2  2er2  In4  2a2 . 

In  the  calculations  of  equilibrium  geometries  vertical  and  adiabatic  excita¬ 
tion  energies  we  used  the  CASSCF  or  state-averaged  CASSCF  methods  [33]. 
In  the  present  work  the  stationary  points  optimization  was  restricted  to  pla¬ 
nar  (D2;,,  D4/i,  C2„)  or  linear  centro-symmetric  structures.  It  has  been  known 
in  previous  calculations  [12]  that  rhombic,  planar  C2,  and  centro-symmetric 
structures  represent  local  minima  on  the  ground  state  PES  of  B4,  B4  and  B4  . 
For  the  B]  linear  structure  the  symmetry  of  the  ground  state  has  been  calcu¬ 
lated  to  be  2 22,  [12]  and  2X„  [31],  whereas  the  present  CASSCF  calculations 
yield  a  2n„  state.  In  the  case  of  the  B4  ground  state  a  2TI„  state  has  been 
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obtained  [12];  we  find  on  the  other  hand  the  lowest  centro-symmetric  linear 
isomer  to  have  the  2n„  symmetry.  The  EOM-CC  method  [34]  was  also  em¬ 
ployed  for  the  B4  excited  singlet  states  and  some  of  the  energy  differences 
were  calculated  by  the  RCCSD(T)  approach  [35].  In  most  calculations  the  basis 
set  used  was  the  augmented  correlation-consistent  polarized  valence  triple  zeta 
(aug-cc-pVTZ)  [36].  The  relative  energies  of  the  isomers,  the  ionization  en¬ 
ergy  and  the  electron  affinity  were  calculated  by  the  RCCSD(T)  method  and 
the  aug-cc-pVQZ  basis  set  [36].  Our  most  reliable  excitation  energies  were 
obtained  by  the  full  valence  CASSCF  method  for  each  state  of  a  given  irre¬ 
ducible  representation  calculated  separately.  For  singlet  states  the  full  valence 
active  space  comprised  about  1.7  x  106  CSF's  for  the  rhombic  and  3.5  x  106 
CSF's  for  the  planar  C2v  structures.  More  information  about  the  computations 
is  specified  in  the  footnotes  of  Tables  1  to  7. 


3.  Results  and  discussion 

The  optimized  D2*  and  C2v  planar  or  centro-symmetric  linear  stationary  points 
on  the  six-dimensional  PES  of  the  electronic  ground  and  excited  states  for  B4, 
B4  and  B4  are  given  in  Tables  1  to  3.  The  geometry  parameters  and  symmetry 
axes  are  defined  in  Fig.  1 .  The  calculated  vertical  and  adiabatic  electronic  exci¬ 
tation  energies  are  compiled  in  Tables  4  to  7.  The  RCCSD(T)  relative  energy 
differences  between  the  optimized  structures  for  the  three  isomers  show  that 
in  all  three  species  the  most  stable  isomers  have  rhombic  structures  followed 
by  planar  C2v  structures  which  lie  2.18  eV  (B4).  1.09  eV  (B4),  1.55  eV  (B4) 
(Tables  1  to  3)  higher  in  energy,  in  relatively  good  agreement  with  previous 
theoretical  findings  (c/.  Ref.  [12]).  The  lowest  dissociation  energy  ( B  1  ( 1 S )  + 
B3  (X  2  A] ))  has  been  determined  experimentally  to  be  2.4  ± 0.6  eV  [6],  though 
the  theoretical  value  was  calculated  several  eV  larger  [17, 29].  The  ground  state 
PES  at  optimized  geometries  of  the  three  isomers  of  B4  lies  below  this  asymp¬ 
totic  limit.  The  rhombus/square  energy  difference  of  0.031  eV  on  the  double 
minimum  potential  of  the  electronic  ground  state  of  B4  leads  to  splittings  of 
vibrational  transitions  [30].  In  B4  (2Ag)  the  barrier  to  the  square  transition 
state  is  calculated  to  be  ten  times  smaller  and  will  not  lead  to  such  splittings 
(Table  2  and  [17]).  In  the  2A„  state  of  B4  this  barrier  is  almost  three  times  larger 
(0.082  eV,  Table  3)  and  thus  the  splittings  will  be  smaller  than  in  rhombic  B4. 
In  the  geometry  region  of  the  square  saddle  point  of  the  B3„  state  the  negative 
ion  is  calculated  to  be  unstable  with  respect  to  autoionization.  The  B2„  state 
correlating  with  the  corresponding  E„  state  autoionizes  already  in  the  geometry 
region  of  the  B3h  minimum. 

In  Bj  the  lowest  2A„  state  lies  only  0.17  eV  below  the  2Blg  state  (Table  6, 
method  e).  Both  states  cross  in  D2*  structures  along  the  ring  angle  coordi¬ 
nate  (Fig.  2a),  and  form  an  avoided  crossing  along  the  b[g  coupling  mode 
one  (Fig.  2b).  The  accuracy  of  the  ah  initio  calculated  energy  differences  and 
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Table  1.  CASSCF  optimized  geometries  for  the  B4 
between  their  minima. 

isomers  and  relative  energy  differences 

B4  ( D2i ,)  rhombus 

B4  (D4;,)  square 

State 

Method 

R 

[bohr] 

a 

[deg.] 

A-£cCSD(T) 

[eV] 

Method 

R  A£'CCSD(T)fl 

[bohr]  [eV] 

X  >AS 

b 

2.909 

75.8 

0 

b 

2.899  0.031 

'B,„ 

c 

3.134 

73.9 

c 

3.266 

'Bif 

c 

3.092 

61.7 

c 

3.033 

'b3„ 

c 

3.007 

69.8 

c 

3.007 

'B3s 

c 

3.058 

60.3 

c 

2.981 

‘B2S 

c 

3.060 

64.7 

c 

2.981 

*b2„ 

c 

3.024 

59.4 

c 

3.007 

1  A„ 

c 

3.300 

90.0 

c 

3.300 

B4  ( C2v ) 

State 

Method 

Ri 

Ri 

e 

^  ^CCSD(T)  ° 

[bohr] 

[bohr] 

[deg.] 

[eV] 

X  ‘A, 

d 

3.227 

3.026 

56.8 

2.182 

e 

3.253 

2.986 

58.4 

‘b2 

f 

3.386 

3.218 

54.4 

IB1 

f 

2.936 

3.100 

58.3 

*a2 

f 

3.207 

3.205 

55.2 

‘A! 

f 

3.111 

3.174 

54.6 

B4  (Do*)  linear 

State 

Method 

r2 

A-£cCSD(T) 

a 

[bohr] 

[bohr] 

[eV] 

X‘X+ 

8 

2.922 

2.890 

2.379 

“  AECCSDm  relative  to  the  X  'A,  state  of  rhombic  B4;  aug-cc-VQZ  basis  set;  E  = 
—  98.91001460  hartree. 
b  CCSD(T)  with  aug-cc-VQZ  basis  set. 

c  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3as,  l-2b3„,  l-2b2„,  l-2b!„,  l-2bls, 

lb2g,  lb3g. 

d  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-6a1;  l-3bj,  l-3b2. 
e  CCSD(T)  with  aug-cc-VTZ  basis  set. 

/  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-7a1;  l-3bl5  l-3b2,  la2. 
g  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3as,  l-3b3„,  l-2b2„,  l-2bln,  lbls,  lb2s. 


geometries  does  not  permit  a  definite  answer  about  the  mixing  of  both  states, 
i.e.  about  the  structure  of  the  vibronic  ground  state,  vibronic  transitions  and 
intensity  borrowing  in  the  UV  absorption  spectra. 
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Table  2.  CASSCF  optimized  geometries  for  the  B4  isomers  and  relative  energy  differ¬ 
ences  between  their  minima. 


b4+  (o2h 

)  rhombus 

B+  (D4„)  square 

State 

Method 

R 

a  A2sRCCSD(T)fl 

Method 

R 

A-ErccSDCT)0 

[bohr] 

[deg.]  [eV] 

[bohr] 

[eV] 

X  2Ag 

b 

2.958 

79.5  0 

b 

2.913 

0.003 

c 

2.951 

83.4 

2Bln 

d 

3.026 

76.1 

d 

3.003 

2Bis 

d 

3.036 

61.6 

d 

2.954 

2b3„ 

d 

3.127 

54.2 

d 

3.030 

2b2„ 

d 

N.C. 

N.C. 

d 

3.030 

2B3g 

d 

3.060 

77.6 

d 

3.155 

2B2g 

d 

3.152 

65.3 

d 

3.155 

2A„ 

d 

3.247 

64.4 

d 

3.159 

B4+  (C2„) 

State 

Method 

Ri 

e 

^  ^RCCSD(T)  ° 

[bohr] 

[bohr] 

[deg.l 

[eV] 

X  2A, 

e 

3.113 

3.087 

56.5 

1.087 

c 

3.086 

3.039 

57.3 

2b2 

f 

3.602 

3.149 

57.1 

2B, 

f 

2.978 

2.986 

121.2 

2A, 

f 

3.081 

3.017 

61.5 

2B, 

f 

3.157 

3.169 

54.1 

2b2 

f 

3.187 

3.196 

56.7 

b4+ 

linear 

State 

Method 

Ri 

Ri 

A  RCCSD(T) 0 

[bohr] 

[bohr] 

[eV] 

x2n„ 

8 

2.895 

2.998 

1.496 

“  AERccsd(T)  relative  to  the  X  2Ag  state  of  rhombic  B4;  aug-cc-VQZ  basis  set;  E  = 
—  98.54990461  hartree. 

b  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3ag,  l-2b3„,l-2b2„,  l-2big,  l-2bi„. 
c  RCCSD(T)  with  aug-cc-VTZ  basis  set. 

d  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3ag,  l-2b3„,  l-2b2„,  l-2blB,  l-2big, 
I  b2,,,  lb3g;  the  optimization  for  the  2B2„  state  did  not  converge. 
e  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-7a1;  l-2bi,  l-3b2. 

/  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-7ai,  l-3b1;  l-3b2,  la2. 
g  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3ag,  l-3b3„,  l-2b2„,  l-2bi„,  lbig,  lb2g. 
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Table  3.  CASSCF  optimized  geometries  for  the  B4  isomers  and  relative  energy  differ¬ 
ences  between  their  minima. 

State 

B4  (Dlh)  rhombus 

B4  ( D4i ,)  square 

Method 

R 

[bohr] 

a 

[deg.] 

^  ^RCCSD(T) 0 

[eV] 

Method 

R  AERcCSD(T)a 

[bohr]  [eV] 

X  2B3„ 

b 

2.992 

63.3 

0 

b 

2.922 

2A, 

c 

3.035 

68.6 

0.050 

c 

3.007  0.082 

B4  (C2v) 

State 

Method 

Ri 

r2 

0  A£RCCSD(T)a 

[bohr] 

[bohr] 

[deg.] 

[eV] 

X  2A, 

d 

3.226 

3.050 

56.8 

1.546 

2B, 

e 

3.097 

3.069 

60.9 

2.014 

B4  linear 

State 

Method 

A, 

Ri 

A£RCCSD(T)a 

[bohr] 

[bohr] 

[eV] 

x2n„ 

/ 

2.919 

3.103 

2.866 

“  A£rccsd(T)  relative  to  the  X  2B3„  state  of  rhombic  B4  ;  aug-cc-VQZ  basis  set;  E  = 
—  98.96948401  hartree. 

b  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-3as,  l-3b3„,  l-2b2„,  l-2b,s,  l-2bln. 

c  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-4a,,  l-2b3„,  l-2b2„,  l-2b,s,  l-2bln. 

d  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-7a1;  l-4bj,  lb2. 

e  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-6aj,  1-4  b,,  l-3b2. 

/  CASSCF  with  aug-cc-VTZ  basis  set;  active  space:  l-4as,  l-3b3„,  l-2b2„,  1  -2b  L„ ,  lblg,  lb2s. 


Similar  coupling  is  found  also  for  the  ground  state  of  B4 .  In  the  adiabatic 
picture  the  2B3li  state  is  calculated  only  0.05  eV  below  the  2Ag  state  (Table  3, 
and  Ref.  [5]).  Also  in  this  case  both  states  cross  in  D2;,  structures  along  the  ring 
angle  coordinate  (Fig.  3a),  and  form  an  avoided  crossing  along  the  b3„  coupling 
mode  (Fig.  3b). 

The  adiabatic  aug-cc-pVQZ/RCCSD(T)  ionization  energy  IE,,  of  B4  has 
been  calculated  to  be  9.799  eV  for  the  rhombic,  8.705  eV  for  the  planar  C2„ 
and  8.916  eV  for  the  linear  structures  using  the  optimized  geometries  (Tables  1 
to  2).  The  first  is  in  excellent  agreement  with  the  experimental  value  of 
9.8  eV  [2-4].  Similar  approach  yields  adiabatic  electron  affinities  of  1.618  eV 
(negative  ion  in  2B3„  state)  and  1.568  eV  (in  2A„  state)  for  the  rhombic  isomer, 
1.786  eV  (in  2B!  state)  and  2.254  eV  (in  2A , )  for  the  C2„  isomer,  and  1.131  eV 
(in  2FI„  state)  for  the  linear  B4  (c/.  also  Ref.  [12]).  The  two  higher  isomers 
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Fig.  1.  Definition  of  the  coordinates  and  symmetry  axis  in  B4  isomers. 


have  well  defined  electronic  wavefunctions  indicating  that  spontaneous  elec¬ 
tron  detachment  does  not  take  place.  The  planar  C2„  isomer  (2Aj)  is  calculated 
to  lie  just  below  the  autodetachment  threshold  of  1.6  eV  [5]  (Table  3).  The  elec¬ 
tronically  excited  states  will  autoionize,  but  the  ground  state  six-dimensional 
PES  has  large  portions  on  which  the  additional  electron  is  bound,  particularly 
along  the  isomerization  paths. 

The  full  valence  CASSCF  calculations  of  vertical  excitation  energies  con¬ 
firm  that  the  ionization  of  B4  leads  to  triplet  states  of  B4  (Table  4).  Our  results 
agree  well  with  the  RCCSD(T)  values  [5],  and  with  the  positions  of  the  first  two 
observed  bands  in  the  photoelectron  spectrum.  This  good  agreement  with  ex¬ 
periment  allows  us  to  assume  that  other  excitation  energies  calculated  with  this 
approach  will  be  of  similar  accuracy. 

The  vibronic  coupling  in  B4  ground  state  mixes  gerade  and  ungerade  wave- 
functions.  In  its  photoelectron  spectrum  only  the  ionizations  to  ungerade  triplet 
states  were  assigned.  In  Table  4  we  compare  the  vertical  excitation  energies 
relative  to  the  singlet  ground  state  of  B4  but  calculated  for  the  optimized  ge¬ 
ometries  of  the  2B3„  and  2A?  states  of  B4 .  Whereas  the  energy  differences 
for  the  first  two  triplet  states  agree  nicely  with  experiment  for  the  geom¬ 
etry  of  the  2B3„  state,  they  disagree  for  the  second  structure.  This  supports 
the  assumption  made  in  [5]  that  the  lowest  part  of  the  PES  for  the  ground 
state  of  the  negative  ion  is  dominated  by  the  2B3„  wavefunction.  The  lowest 
triplet  states  3B2„  and  3B3„  correlate  with  the  3E„  state  in  D4/,  symmetry;  the 
observed  band  corresponds  to  a  Jahn-Teller  system.  The  experimental  spec¬ 
trum  at  higher  energies  is  more  difficult  to  assign  because  several  triplets 
and  singlets  are  calculated  to  lie  close  to  the  observed  band  at  about  2.4  eV 
(Table  4). 

In  Table  5  the  vertical  and  adiabatic  transition  energies  in  rhombic  and 
square  geometries  for  the  singlet  states  of  B4  are  compared.  Only  one  elec- 
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Table  4.  Comparison  of  the  experimental  and  calculated  energy  differences  in  B4  at  the 
optimized  geometries  of  the  2B,„  and  2AS  states  of  rhombic  B4 . 


'~p  a 

1  V 

[eV] 

Ty 

[eV] 

Tve 

[eV] 

Experiment 

[eV] 

X  ‘As 

0 

0 

0 

0 

a  3B2„ 

1.29 

2.18 

1.29 

1.09 

b  3B3„ 

1.43 

1.95 

1.39 

1.42 

c3Bls 

1.64 

1.84 

1.74d 

A‘Bls 

2.03 

B  'B2„ 

2.13 

d3As 

2.55 

1.75 

1 .60" 

e  3B2? 

2.60 

2.80 

2.50 

C  ‘B2g 

2.84 

f3B,„ 

2.85 

2.37 

2. 1 8rf 

2.40 

D  ‘Bj„ 

3.36 

“  Full-valence  CASSCF;  aug-cc-pVTZ  basis  set;  each  state  calculated  separately  at  the 
optimized  geometry  of  the  X  2B3„  state  of  B4:  R  =  3.015  bohr,  a  =  63.6  degrees  [5]; 
E  =  —98.69773452  hartree. 

b  Full-valence  CASSCF;  aug-cc-pVTZ  basis  set;  each  state  calculated  separately  at  the 
optimized  geometry  of  the  X  2AS  state  of  B4:  R  =  3.064  bohr,  a  =  73.8  degrees  [5]; 
E  =  —  98.69847200  hartree. 

1  RCCSD(T);  each  state  calculated  at  the  optimized  geometry  of  the  X  2B3„  state  of  B4 : 
R  —  3.015  bohr,  a  =  63.6  degrees  [5], 

d  RCCSD(T);  each  state  calculated  at  the  optimized  geometry  of  the  X  2Ag  state  of  B4 : 
R  =  3.064  bohr,  a  =  73.8  degrees  [5], 


TableS.  Vertical  (Tv)  and  adiabatic  ( Te )  excitation  energies  and  transition  moments  (fi,,) 
between  singlet  states  of  rhombic  and  square  B4. 


Rhombus  Square 


'~p  Cl 

1  v,av 

[eV] 

E-e 

[debye] 

Ty,  EOM^ 

[eV] 

V 

[eV] 

rp  d 

[eV] 

'~p  a 

-*  v,av 

[eV] 

E-e 

[debye] 

rp  b 

1  v,EOM 

[eV] 

V 

[eV] 

[eV] 

X 

‘A, 

0 

0 

0 

0 

0 

0 

0 

0 

A 

‘Big 

2.91 

3.12 

2.94 

2.12 

2.98 

3.35 

3.10 

2.87 

B 

'b2„ 

3.44 

0.6889  (v) 

3.71 

3.52 

2.20 

4.57 

1.6642  (y) 

4.63 

4.56 

4.19 

C 

1  B,„ 

3.56 

0.2412  (z) 

3.49 

3.61 

2.89 

3.40 

0.0000  (z) 

3.43 

3.46 

3.30 

D 

'B2g 

3.74 

3.74 

3.77 

2.94 

5.19 

4.48 

5.30 

4.68 

E 

‘b3„ 

3.99 

2.3070  (x) 

3.85 

4.12 

3.82 

4.57 

1.6642  (x) 

4.63 

4.56 

4.19 

F 

'B3g 

4.66 

4.55 

4.58 

3.90 

5.19 

4.48 

5.30 

4.68 

“  Full-valence  SA-CASSCF  with  aug-cc-pVTZ  basis  set  at  the  optimized  geometry  of 
X  1  A,,  (Table  1);  one  state  per  irreducible  representation  calculated. 
b  At  the  optimized  geometry  of  X  1  Ag  (Table  1);  aug-cc-pVQZ  basis  set. 
c  Full-valence  CASSCF  with  aug-cc-pVTZ  basis  set  at  the  optimized  geometry  of  X  1  A,, 
(Table  1);  each  state  calculated  separately. 

d  Full-valence  CASSCF  with  aug-cc-pVTZ  basis  set  at  optimized  geometry  (Table  1);  each 
state  calculated  separately. 
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Table  6.  Vertical  (Tv)  and  adiabatic  ( Te )  excitation  energies  and  transition  moments  (/x,) 
between  doublet  and  quartet  states  of  rhombic  and  square  B4  . 


Rhombus  Square 


y  a 

*  v,av 

[eV] 

[debye] 

T‘ 

[eV] 

y  d 

[eV] 

rp  a 

v,av 

[eV] 

/*.* 

[debye] 

Tvc 

[eV] 

[eV] 

X2Ag 

0 

0 

A  2Blg 

0.58 

0.70 

0.17 

0.73 

0.78 

0.74 

B  2B,„ 

0.63 

0.9796  (z) 

0.70 

0.62 

0.68 

1.0374  (z) 

0.75 

0.60 

a4Bls 

1.99 

2.13 

b4B„, 

2.35 

2.35 

C2AS 

2.50 

c4B2„ 

2.72 

1.4404  (x) 

3.94 

1.1402  (x) 

D  2B1„ 

2.93 

0.2431  (z) 

E2Bis 

2.94 

F  2B3„ 

3.03 

1.7305  (x) 

3.04 

1.64 

3.69 

0.7120  (x) 

3.52 

3.35 

G  2B2„ 

3.17 

1.3820  (y) 

3.69 

0.7120  (y) 

3.52 

3.35 

d  4B2g 

3.21 

4.39 

e  4A„ 

3.28 

0.8564  (z) 

3.74 

0.9007  (z) 

H  2B2g 

3.63 

3.71 

2.46 

4.96 

4.79 

3.75 

f4B3g 

3.75 

4.39 

I  2b2„ 

3.81 

0.6198  (y) 

J  2A„ 

4.07 

4.15 

2.93 

4.68 

4.44 

3.48 

g  4b3„ 

4.13 

0.1174  (y) 

3.94 

1.1402  (y) 

K  2A„ 

4.17 

L  2B3„ 

4.28 

0.8819  (x) 

M  2B3g 

4.31 

4.44 

3.86 

4.96 

4.79 

3.75 

“  Full-valence  SA-CASSCF  with  aug-cc-pVTZ  basis  set  at  the  optimized  geometry  of 
X  2Ag  (Table  2);  two  states  per  irreducible  representation  calculated. 
b  For  the  quartet  states  the  transition  moments  are  given  with  respect  to  the  lowest  4B,g. 
c  Full-valence  CASSCF  with  aug-cc-pVTZ  at  the  optimized  geometry  of  X  2Ag  (Table  2); 
each  state  calculated  separately. 

d  Full-valence  CASSCF  with  aug-cc-pVTZ  basis  set  at  optimized  geometry  (Table  2);  each 
state  calculated  separately. 


tronic  state  in  a  given  irreducible  representation  has  been  considered.  The 
state  averaged  calculations  needed  for  the  transition  moments  yield  tran¬ 
sition  energies  which  deviate  from  those  calculated  for  each  state  sepa¬ 
rately  by  less  than  about  0.1  eV.  The  CC-EOM  vertical  transition  energies 
differ  from  the  full  valence  CASSCF  results  by  up  to  0.2  eV.  The  ab¬ 
sorption  spectrum  of  B4  will  have  strong  bands  (cf.  transition  moments 
in  Table  5)  to  the  'B2,,  state  starting  around  2.20  eV  and  the  'B3„  state 
around  3.82  eV.  These  states  form  a  Jahn-Teller  pair,  the  D4/l  E„  structure 
is  calculated  to  lie  2  eV  and  1 .45  eV  above  the  equilibrium  energies  of  the 
two  rhombic  components.  Though  the  transitions  to  the  Jahn-Teller  pair 
in  gerade  symmetry  are  symmetry  forbidden,  they  could  gain  intensity  by 
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Table  7.  Vertical  (Tv)  and  adiabatic  ( Te )  excitation  energies  and  transition  moments  (fie) 
between  singlet  states  of  C2„  B4  and  between  doublet  states  of  C2„  B4 . 


B4  B+ 


rp  a 
v,av 

[eV] 

/A 

[debye] 

rp  b 

1 1; 

[eV] 

t; 

[eV] 

rp  a 
-*  v,av 

[eV] 

/A 

[debye] 

rp  b 

1 1; 

[eV] 

Tec 

[eV] 

X 

‘A, 

0 

0 

0 

X2A,  0 

A 

‘A, 

1.27 

1.6832  (x) 

1.35 

0.41 

A  2B,  1.33 

0.3485  (z) 

1.37 

1.15 

B 

‘B, 

1.90 

0.9179  (y) 

1.76 

1.50 

B  2 A[  1.47 

0.2443  (jc) 

1.52 

1.32 

C 

'b2 

2.37 

0.2348  (z) 

2.21 

2.01 

C  2B,  1.71 

0.4057  (y) 

1.76 

1.06 

D 

'b2 

2.43 

1.2223  (z) 

N.C. 

D  2Ai  1.94 

E 

■a2 

2.50 

2.53 

2.35 

E  2B,  2.04 

0.6730  (z) 

2.03 

1.93 

F 

‘B, 

2.58 

0.3910  (>•) 

2.61 

F2B,  2.51 

0.5600  (>•) 

2.58 

2.51 

G 

*a2 

3.10 

3.75 

G  2Ai  2.73 

“  Full-valence  SA-CASSCF  with  aug-cc-pVTZ  basis  set  at  the  optimized  geometry  of 
X  'A,  (Table  3);  two  states  per  irreducible  representation  calculated. 
b  Full-valence  CASSCF  with  aug-cc-pVTZ  basis  set  at  the  optimized  geometry  of  X  'A, 
(Table  3);  each  state  calculated  separately. 

c  Full-valence  CASSCF  with  aug-cc-pVTZ  basis  set  at  optimized  geometry  (Table  3);  each 
state  calculated  separately. 


vibronic  coupling  via  the  b3li  mode  with  the  close  lying  ungerade  states 
(Table  5). 

In  B|  two  excited  doublet  states  lie  lower  than  the  first  quartet  state 
(Table  6).  As  discussed,  the  2B  „  state  is  calculated  only  0.17  eV  above  X  2A,, 
and  both  are  vibronically  coupled.  Due  to  possible  intensity  borrowing  it  is  dif¬ 
ficult  to  calculate  reliable  radiative  transition  probabilities  for  the  absorption 
spectra.  Nevertheless  in  the  doublet  states  the  most  intense  transitions  should 
again  be  those  to  the  Jahn-Teller  states  2B3„  and  2B2„  for  which  the  vertical 
energy  difference  is  calculated  to  be  about  3.1  eV.  The  computations  of  the  sta¬ 
tionary  point  in  D2/i  symmetry  for  the  second  2B2„  component  (Table  2)  did 
not  converge  because  of  the  perturbation  by  other  electronic  states.  Also  the 
vertical  energy  differences  for  the  quartet  states  are  given  in  Table  VI.  The 
transition  moments  between  these  states  are  given  relative  to  the  lowest  4B  ? 
state. 

The  C2„  minima  are  calculated  higher  than  those  of  the  rhombic  isomers 
(Tables  1  and  2).  In  Table  7  the  vertical  and  adiabatic  excitation  energies  for 
these  isomers  are  given.  Providing  that  both  rhombic  and  C2„  isomers  co¬ 
exist  in  the  sample,  the  UV  spectra  will  show  transitions  for  both  species. 
In  B4  the  most  intense  transitions  are  calculated  to  the  'A,,  1 B ,  and  1 B2 
states  (Table  7),  which  lie  much  lower  that  the  first  allowed  singlet  transition 
in  the  rhombic  isomer.  This  pattern  should  allow  the  isomers  to  be  distin¬ 
guished  from  their  electronic  spectra.  In  the  case  of  B^  the  comparison  of 
energy  differences  in  Tables  6  and  7  shows  that  the  pattern  of  electronic  tran- 
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B 


+ 

4 


b)  AQBlj/bohr 

Fig.  2.  One  dimensional  cut  in  D2;,  symmetry  along  the  ring  angle  coordinate  and  R,,  = 
2.958  bohr  (cf.  Table  2)  of  the  2Ag  and  2Blg  states  of  Bj.  Lower  (broken  line)  and 
upper  (full  line)  potentials  correspond  to  the  2Ag  and  2Blg  states  respectively.  In  the 
SA-CASSCF  calculations  (cc-pVTZ,  active  space  3-5ag,  2-3b3„,  2-3b2„,  lblg,  l-2blu,  lb2g, 
lb3g)  both  states  were  averaged  with  equal  weights.  In  lower  part  of  the  figure  the  avoided 
crossing  along  the  B  lg  coupling  mode  for  R'f  =  2.958  bohr  and  apex  angle  of  o'"1  = 
68  degrees  is  shown. 


sitions  will  make  such  an  identification  more  difficult,  though  also  in  this 
case  bands  below  about  1.5  eV  will  advocate  the  presence  of  the  second  iso¬ 
mer. 
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154- 


Fig.  3.  One  dimensional  cut  in  D2/,  symmetry  along  the  ring  angle  coordinate  and  Re  = 
2.992  bohr  (cf.  Table  2)  of  the  2Ag  and  2B3„  states  of  B4 .  Lower  (full  line)  and  upper  (bro¬ 
ken  line)  potentials  correspond  to  the  2B3„  and  2AS  states  respectively.  In  the  SA-CASSCF 
calculations  (cc-pVTZ,  active  space  3-6ag,  2-3b3„,  2-3b2„,  lbls,  lbi„,  lb2s,  lb3s)  both 
states  were  averaged  with  equal  weights.  In  lower  part  of  the  figure  the  avoided  crossing 
along  the  B3„  coupling  mode  for  R'f  =  2.992  bohr  and  apex  angle  of  aref  =  67  degrees  is 
shown. 


4.  Summary  and  conclusions 

So  far  only  the  photodetachment  spectrum  of  rhombic  [5]  has  provided  ex¬ 
perimental  information  on  the  singlet  ground  and  excited  triplet  states  of  B4  for 
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a  comparison  with  theoretical  results.  The  calculated  full  valence  CASSCF  ex¬ 
citation  energies  for  the  observed  two  lowest  triplet  states  agree  to  within  about 
0.1  eV  with  the  experimental  band  positions.  A  similar  computational  approach 
was  employed  to  calculate  vertical  and  adiabatic  excitation  energies  up  to  about 
4.5  eV  for  B4  and  Bj ,  including  the  electronic  transition  moments,  which  could 
aid  the  search  for  boron  tetramers  by  electronic  spectroscopy.  In  the  rhom¬ 
bic  ground  state  of  the  charged  clusters  the  Born-Oppenheinrer  approximation 
fails;  in  B4  the  2A?  and  2Bls  states  are  coupled  by  the  b|„  vibration,  in  B4  the 
2AS  and  2B3„  states  by  the  b3„  mode. 
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